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INTRODUCTION 


2.  EARLIER  WORK  ON  CYCLONE  FREQUENCIES 


Damaging  wavas  and  storm  surges 
along  tha  U.S.  Atlantic  Coast  ara  largely  due  to 
extratropical  cyclones.  Along  moat  of  this  coast, 
shorelines  are  receding  (U.S.  Army  Corps  of  Engi¬ 
neers,  1971;  Dolan  at  al. ,  1979).  The  erosion 
trend  is  attributed  to  tha  currant  rise  in  saa 
level  (Bruun,  1962) ,  lower  average  cyclone  central 
pressure  (Mather,  Adams,  and  Yoshioka,  1964),  re¬ 
duced  coastal  sand  supply  (Hoyt,  1967),  human 
activities  (Dolan,  Godfrey,  and  Odisa ,  1973),  and 
to  secular  variation  in  cyclone  frequency,  magni¬ 
tude  and  duration  (Haydan,  1975).  Mora  racantly, 
changes  in  the  tracks  of  cyclones  along  the  At¬ 
lantic  coast  have  been  reported  (Resio  and  Hayden, 
1975;  Dickson  and  Namias,  1976).  Changes  in  the 
track  of  cyclones  relative  to  the  coast  give  rise 
to  variations  in  breaker  heights  and  storm  surge 
heights  at  the  coast  (Resio  and  Haydan,  1975). 

The  Resio  and  Haydan  study  is 
limited  by  the  small  reach  of  coast  studied  and 
the  Dickson  and  Namias  study  by  the  short  time 
span  analysed  (1948-1975).  In  this  report,  spa¬ 
tial  and  temporal  variations  in  cyclone  frequenciea 
are  examined  for  the  94  years  between  1885  and 
1978.  The  study  area  is  Eastern  North  America 
and  the  western  North  Atlantic  east  of  the  100th 


Figure  l.  Chart  of  the  study  area.  The  rectan¬ 
gular  ineet  is  2.5*  latitude  by  5.0 •  longitude. 
There  are  74  such  rectangular  grid  sells  in  the 
study  area.  The  arrou  represents  a  stom  track 
passing  through  the  grid  sell  shoun. 


Petterssen  (19S0)  mapped  the  geograph¬ 
ical  frequencies  of  cyclones  in  the  Northern 
Hemisphere  using  daily  sea-level  Historical  Weather 
Maps  for  the  period  1899-1939.  Petterseen's  map 
shows  a  frequency  maximum  along  the  eaet  coest  of 
the  U.S.  as  far  south  as  Florida.  Klein  (1957) 
attributes  the  east  coast  cyclone  frequency  maxi¬ 
mum  to  an  "Atlantic  seaboard"  cyclone  track. 
Cyclogenesis  associated  with  this  track  occurs 
along  the  tier  of  Gulf  Coast  states  or  along  the 
Virginia  Capes.  Reed  (1960)  indicates  in  his  maps 
of  percentage  frequency  of  fronts  that  the  front 
frequency  maxi mum  also  parallels  the  east  coast 
margin  but  at  some  distance  offshore  during  tha 
years  1952-1956.  Miller  (1946)  analyzed  cyclogenesis 
for  the  mid-Atlantic  coastal  region  and  found  two 
distinct  types.  His  Type  &  form  largely  to  the 
north  of  Cape  Hatteras,  North  Carolina ,  and  Type 
B  to  the  south  of  the  Cape .  while  seasonal  vari¬ 
ations  in  cyclone  frequency,  front  frequency  and 
track  locations  are  discussed  in  some  detail, 
temporal  variations  ware  not  studied. 

More  recent  studies  of  cyclones  along 
the  Atlantic  coast  have  focused  on  the  generation 
of  potentially  damaging  waves  and  surges.  Bosser- 
man  and  Dolan  (1968)  show  an  increase  in  the  fre¬ 
quency  of  storms  that  generate  deep-water  waves 
of  1.3  m  or  greater  from  1942  to  1967.  The  Boaser- 
man  and  Dolan  study  was  updated  by  Haydan  (1975) .  • 

Haydan  concluded  that  the  number,  severity  and 
duration  of  Atlantic  coast  storms  had  increased 
between  1942  and  1974.  Increases  in  cyclone  sever¬ 
ity  are  also  reported  by  Mather  et  al. ,  (1964)  and 
they  attribute  tha  changes  to  decreased  central 
pressures  between  the  1920's  and  1960's.  Resio 
and  Hayden  (1975)  found  that  cyclones  tracking 
along  the  U.S.  Gulf  Coast  and  the  genesis  of  new 
cyclones  along  the  mid-Atlantic  coast  had  increased 
from  the  1920's  to  the  196Q's.  The  pre-1920's  period 
more  closely  resembled  tha  I960' s  than  the  years 
following  tha  1920' s.  Resio  and  Hayden  also  zhowed 
that  the  increased  severity  of  Atlantic  coastal 
storms  resulted  from  a  seaward  displacement  of 
mean  storm  tracks.  They  attributed  this  Change  to 
increased  blocking  in  the  high  latitudes.  Dickson 
and  Namias  (1976)  identified  blocking  in  the  Green¬ 
land  area  and  cooler  than  normal  t«>pe returns  in 
the  U.S.  Southeast  as  the  cause  of  increased  fre¬ 
quencies  of  cyclones  off  the  Atlantic  coast  between 
the  1940 ' s ,  and  the  1950's  and  1960's. 
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3. 


THE  DATA 


Cyclone  frequencies  for  *ach  year 
(1385-19781  vara  tabulated  for  the  74,  2 .  5  "-lati¬ 
tude  by  5.0 "-longitude  grid  call*  comprising  tha 
study  araa  shown  in  Figure  1.  From  monthly 
charts  of  tha  “Tracks  of  tha  Cantaxs  of  Cyclonas 
at  Saa  Laval"  publishad  by  tha  Monthly  Weather 
Review  and  in  racant  yaars  by  tha  Mariners 
Weather  Log,  annual  totals  of  cyclonas  passing 
through  each  grid  call  wara  recorded.  Multiple 
entries  of  a  given  star*  in  a  grid  call  ware 
ignored.  Tha  resulting  data  matrix  has  tha 
dimensions  of  74  grid-call- variable*  and  94 
annual  cases. 

4.  LONG-TERM  MEANS  AND  VARIANCES 

figure  2  shows  tha  mean  annual 
frequency  of  cyclones  for  tha  period  1885-1978. 

In  general  tha  frequency  of  cyclones  increases 
with  latitude  and  with  proximity  to  tha  O.S. 
east  coast.  Tha  axis  of  tha  east  coast  fre¬ 
quency  maxima  is  centered  over  tha  coast  from 
Move  Scotia  to  Georgia.  Tha  location  of  this 
frequency  maxima  coincides  with  tha  baroclinic 
zona  separating  tha  cold  mainland  from  tha  warmer 
ocean  waters. 

Tha  standard  deviation  of  annual 
cyclone  frequencies  is  presented  in  Figura  3. 

Tha  dominant  feature  shown  by  tha  chart  is  tha 
pronounced  frequency  maximum  offshora.  Tha  axis 
of  this  maximum  is  seaward  of  the  maximum  in  the 
mean  field  (Fig.  2)  and  apparently  rotated  clock¬ 
wise  several  degrees.  -Tha  location  of  tha  stan¬ 
dard  deviation  maxlsnsa  more  closely  approximates 
tha  locus  of  tha  Gulf  Stream  than  the  continental 
margin.  Tha  standard  deviations  landward  of  the 
margin  of  the  cast  coast  are  rather  modest  sug¬ 
gesting  that  tha  axis  of  yearly  storm  occurrences 
may  expand  eastward  but  in  general  not  over  tha 
mainland. 


Figure  2.  Hoax  annual  oyalona  fraquanoy  for  tha 
yaara  1885-1978. 


Figura  3.  Standard  daviationa  of  tha  moan  annual 
oyalona  fraquaney  for  tha  yaara  1885-1978. 


S.  PRINCIPAL  COMPONENTS  ANALYSIS  OP 

ANNOAL  FREQUENCY  DATA 


The  Anely 


Principal  coapanent  analysis  has 
successfully  resolved  tha  variance  structure  in 
multi variant,  geophysical  data  (Xutxbach,  1967) 
Frltts  at  el.,  1971;  Rasio  and  Hayden,  1975),  and 
in  terms  of  least  square  errors,  this  type  of 
analysis  provides  a  method  for  determining  pat¬ 
terns  in  large  data  fields  (Lorenz,  1956; 

Gilman,  1957;  Kutzbach,  1967).  Tha  abjective  of 
tha  analysis  ia  to  isolate  characteristic,  recur¬ 
rent,  and  independent  modes  of  covariance  among 
variables  into  a  new  sat  of  independent  variables. 
Basically,  the  snelymle  transforms  a  sat  of  inter- 
correlated  variables  into  a  new  coordinate  system 
in  which  the  exes  ere  linear  eaafeinatione  of  the 
original  variataa  end  arm  mutually  orthogonal. 


To  prevent  thoee  grid  cells  with 
high  mean  cyclone  frequencies  (high  latitudes) 
from  dominating  tha  total  variance  and  conse¬ 
quently  from  dominating  the  eigenvector  forme, 
the  correlation  matrix  was  used  rather  than  the 
covariance  matrix.  Tha  procedure e  far  calculation 
of  tha  eigenvectors  follow  those  of  Kutxbach  (1967) 
and  Vincent  et  el.  (1976) . 


Principal  component  analysis  pro¬ 
vides  t  description  of  the  major  modes  of  vari¬ 
ability  la  tha  data  set.  Typically,  each  component 
ia  identified  with  acme  property  of  the  data  field. 
The  analysis  also  provides  an  index  which  measures 
the  importance  of  each  component  within  each  year. 
Finally,  the  analysis  provides  an  estimate  of  the 
total  percent  of  variance  in  the  data  set  which 
can  be  explained  on  the  basis  of  each  component. 


The  jU 


The  percentage  of  variance  and  the 
cumulative  percentage  of  variance  explained  by  the 
first  four  eigenvectors  is  given  in  Table  1.  As 
indicated  ia  the  table  the  first  four  eigenvectors 
account  for  57.4%  of  tha  total  variance.  The  prob¬ 
lem  has  thus  been  reduced  from  e  74  variable 
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problem  in  which  the  variables  are  intercorre- 
lated,  to  a  4  variable  problem  in  which  each  new 
variable  ia  orthogonal  and  thus  statistically 
independent.  Higher  order  eigenvectors  explained 
smaller  fractions  of  the  total  variance  and  are 
not  dealt  with  here. 


Eigenvector 

Number 


%  Variance 
Explained 


Cumulative  % 
Variance  Explained 


Table  1.  The  percentage  of  the  total  variance 
and  oumlative  percentage  of  the  total  variance 
of  annual  cyclone  frequencies  associated  with 
eigenvectors  1-4.  The  values  in  column  two  are 
equal  to  the  corresponding  eigenvalue  divided 
by  the  sum  of  all  74  eigenvalues. 

The  eigenvectors  corresponding  to 
the  eigenvalues  1  through  4  are  mapped  in  Figures 
4,  5,  6,  and  7.  The  first  eigenvector  {Fig.  4) 
indicates  that  the  dominant  mode  of  cyclone  fre¬ 
quency  variation  about  the  mean  is,  with  positive 
weightings,  an  abundance  of  cyclones  over  marine 
areas  at  the  expense  of  cyclone  frequencies  over 
the  continent.  Years  with  negative  weightings 
would  be  characterised  by  few  cyclones  over  ma¬ 
rine  areas  and  above  average  frequencies  over 
the  continent.  The  corridor  of  positive  values 
over  marine  areas  shown  in  Figure  4  corresponds 
in  location  to  the  zone  of  increased  cyclone  fre¬ 
quencies  found  by  DidVson  and  Namias  (1976)  and 
the  east  and  south  displaced  storm  track  type 
identified  by  Hasio  and  Hayden  (1975) . 


Figure  4.  The  first  eigenvector  of  annual 
cyclone  frequencies. 

The  second  eigenvector  (Fig.  5) 
has  positive  values  over  the  entire  field  and  as 
such  indicates  for  positively  weighted  years,  a 
general  increase  in  cyclone  numbers  and  in  nega¬ 
tively  weighted  years,  fewer  than  average  cy¬ 
clones.  The  dominant  feature  of  the  map  of 
eigenvector  no.  2  is  the  mar  lam  centered  over 
the  coast  of  the  mid-Atlantic.  This  pattern  is 
rather  similar  to  the  chart  of  cyclogenesis  fre¬ 
quency  published  by  Pettersaen  (1941). 


This  suggests  that  the  second  principal  compo¬ 
nent  might  best  be  termed  a  coastal  cyelogenesis 
function. 


Figure  S.  The  second  eigenvector  of  annual  cy¬ 
clone  frequencies. 

The  third  eigenvector  (Fig.  6)  con 
trasts  storms  tracking  eastward  from  the  Great 
Plains  up  the  Ohio  Valley  and  exiting  the  coast 
between  Cape  Hattaras  and  Cape  Cod,  with  frequent 
cyclones  in  the  Gulf  Coast  region. 


Figure  8,  The  third  eigenvector  of  cetnual  cy¬ 
clone  frequencies. 

The  fourth  eigenvector  (Fig.  7)  contrasts  the  fre 
quant  cyclone  occurrence  in  the  Great  Lakes-Saint 
Lawrence  River  region  with  cyclones  out  of  the 
Great  Plains.  In  years  with  positive  weightings, 
storsM  are  more  frequent  than  average  in  southern 
Canada.  Negatively  weighted  years  shew  an  in¬ 
crease  in  storm  frequency  in  the  Colorado  storm 
track  area. 

The  first  four  eigenvectors  of 
annual  cyclone  frequencies  constitute  four  new 
orthogonal  axes  which  account  for  nearly  60%  of 
the  variance  in  the  original  data.  In  general, 
weightings  on  these  four  vectors  for  the  94  indl- 
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The  moa ninga  of  each  of  the  four  new  axaa  ara 
tuanariZBd  in  Figure  3. 


Figure  7.  Tha  fourth  eigenvector  of  cornual 
ay clona  fraquanciaa . 


4  NEW  AXES  OF  CYCLONE  FREQUENCY 
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Figure  9.  Tim  variation  of  the  annual  weightings 
of  the  fivet  eigenvector  (E^i . 

Annual  weightings  of  tha  tacand 
•igen vector  (Fig.  10)  also  axhibit  a  can tury- Lang 
scala  ot  variation  ffrcxn  nagativa  valuaa  at  tha 
baginning  of  tha  time  sariaa  to  positive  valuaa 
aftar  1910  and  apparently  hack  to  nagativa  valuaa 
around  1960.  Xf  tha  interpretation  of  aiganvaetor 
2  as  ropra oanting  coastal  cydoganaais  is  correct, 
than  thara  should  hava  baan  an  incraaaa  in  ooaatal 
cydoganaais  during  tha  first  fiva  dacadas  of  tha 
20th  can tury .  This  obaarvation  is  supported  by 
Basic  and  Haydan  (1979)  who  show  an  incraaaa  in 
cydoganaais  along  tha  mid-Atlantic  coast. 
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Figure  8.  SteransaXvon  of  the  locations  of 
cyclone  numbers  in  exeeee  of  the  mean  for  each 
of  the  first  four  eigenvectors.  The  -  sign 
aid  the  *  sign  associated  with  ayologeneeis  for 
axis  no.  2  indicates  below  and  above  normal 
cyclone  murders  in  the  region  of  coastal 
ayalogenesis  along  the  U.S.  mid-Atlantia  Coast. 

5.3  Tha  Tima  Sorias 

Tha  tiaa  sariaa  of  aiganvaetor 
weightings  tor  sach  yaar  (1885-1978)  for  aaeh  of 
tha  first  four  eigenvectors  is  show  in  Figures 
9,  10,  11,  and  12.  Tha  annual  weightings  of  tha 
first  aiganvaetor  show  a  broad  gsasi-oinusoidal 
variation  of  period  slightly  greater  than  100 
years.  This  variation  thus  suggests  that  cyclone 
frequencies  dec lined  ever  marina  areas  and  la- 
crossed  over  the  continent  batsman  13 IS  and  1925. 
Sines  1925  eyelone  frequencies  over  sarins  areas 
hava  inezeased  and  frequencies  over  the  continent 
hava  daeliaad.  With  tha  11  si  rad  tiaa  so  rise 
length  it  is  not  possible  to  establish  whether 
tha  variation  ia  indeed  periodic. 


Figure  10.  Tim  variation  of  the  actual  weightings 
of  the  second  eigenvector  (E^J. 

Tha  time  sariaa  of  tha  third  and 
fourth  eigenvectors  (Fig.  11  and  Fig.  12)  show 
excursions  from  positive  to  negative  weightings 
over  shorter  periods  of  years.  Tha  third  eigen¬ 
vector  cuss  from  nagativa  valuaa  to  positive 
values  between  1900  and  1929,  1929  ad  1990,  and 
between  1999  and  1970.  Spootxal  analysis  shows 
sons  concentration  of  spnrtral  power  at  abeat  23 
years.  The  fourth  aiganvaetor  shows  a  bread 
dentins  from  positive  veleas  is  tha  early  past  of 
tha  record  to  lew  values  is  tha  lose  1940's, 
tinea  tha  lata  1940's  weightings  hava  again  in- 
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Figure  11.  Time  variation  of  the  annual 
’-•righting 8  of  the  third,  eigenvector  (S^S . 
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The  variance  in  cyclone  frequency 
is  maximum  off  the  U.S.  Atlantic  Coast.  Much  of 
this  variance  is  due  to  the  eastward  displacement 
of  storm  tracks  that  parallel  the  coast  and  to 
cyclogenesis  along  the  mid-Atlantic  portion  of 
the  east  coast.  Both  of  these  two  contributions 
to  the  total  variance  exhibit  century-long  secular 
variation.  These  changes  are  undoubtedly  due  to 
changes  in  the  intensity  of  the  east  coast  baro- 
clinic  zone  and  to  shifts  in  the  North  American 
long-wave  location  associated  with  blocking  in 
the  high  latitudes  as  suggested  by  Resio  and  Hay¬ 
den  (1975)  and  Dickson  and  N ami  as  (1976) . 

It  thus  appears  that  the  increased 
storminess  of  the  U.S.  east  coast  sines  the  early 
1940 's  (Mather  et  al. ,  1964;  Bosserman  and  Dolan, 
1968;  Hayden,  1975;  Rssio  and  Haydan,  1975; 

Dickson  and  Namias,  1976)  is  part  of  a  secular 
variation  of  longer  time  scale. 
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ABSTRACT 


Winter  and  summer  cyclone  frequencies  for  eastern 
North  America  and  the  western  North  Atlantic  were 
tabulated  for  2.5°  latitude  by  5°  longitude  grid 
cells  for  the  years  1885-1979.  Correlation  matrix 
eigenvectors  were  calculated  for  matrices  of  both 
the  winter  and  summer  data.  The  first  four  eigen¬ 
vectors  of  each  matrix  are  highly  similar  in  form  and 
the  first,  second  and  fourth  eigenvectors  are  highly 
correlated  in  the  time  domain.  These  correlations 
permit  a  season-in-advance  estimation  of  cyclone 
frequencies,  '’’he  multivariate  prediction  model  developed 
is  called  the  University  of  Virginia  Climate  Prediction 
Model,  Forecast  skill  in  both  winter  and  summer 
averages  66%  relative  to  the  means  as  forecasts. 

Forecast  failure  (negative  skill)  occurs  in  about  8  percent 
of  the  winter  forecasts  and  11  percent  of  the  sumner  forecasts. 
Positive  skill  relative  to  chance  and  persistence 
is  also  demonstrated. 
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1.  Introduction 

In  1914,  C.J.  Kullmer  proposed  that  cyclone  frequency 
patterns  for  North  America  were  in  part  determined  by  the 
motions  of  sunspots  on  the  solar  surface  and  were  thus 
predictable  to  a  degree.  Kullmer' s  subsequent  work  (1933, 
1943)  further  substantiated  his  claims.  Ellsworth  Huntington 
(1914a,  1914b)  publicized  Kullmer's  ideas  as  the  "Solar 
Cyclonic  Hypothesis."  The  Solar  Cyclonic  Hypothesis  has 
since  fallen  into  obscurity  but  the  goal  of  long-range 
climate  prediction  has  been  codified  in  the  National 
Climate  Program.  In  this  report  a  method  for  season- in¬ 
advance  cyclone  frequency  prediction  is  detailed.  The 
multivariate  statistical  model  developed  is  referred  to  here 
as  the  University  of  Virginia  Climate  Prediction  Model. 

2.  Earlier  Work 

a.  Cyclone  Frequencies 

The  movements  of  the  centers  of  cyclones  have  been 
charted  for  North  America  since  1873.  Systematic  study 
of  the  climatology  of  cyclones  began  in  the  1890s  with 
the  work  of  Bigelow  (1897)  on  "storms,  storm  tracks  and 
weather  forecasting."  Shortly  after  the  turn  of  the  century 
two  groups  began  detailed  studies  of  cyclones  and  cyclone 
movements.  Bowie  (1906)  and  Bowie  and  Weightman  (1914) 
studied  the  frequency  of  cyclones  geographically  and 
classified  the  cyclones  by  origin  and  movement  in  order 
to  improve  weather  forecasting.  The  other  group,  Kullmer 


(1914)  and  Huntington  (1914a,  19140,  studied  covariations  in 
the  record  of  sunspots  and  cyclone  frequencies  in  order 
to  achieve  long-range  cyclone  frequency  prediction. 

Bowie  and  Weightman's  work  has  stood  the  test  of  time. 
Shortly  thereafter  Bjerknes  (1918)  deduced  the  three- 
dimensional  structure  of  moving  cyclones  and  with  Solberg 
(Bjerknes  and  Solberg,  1922)  published  a  paper  on  the  life  cycle  of 
cyclones  and  the  theory  of  the  polar  front.  In  contrast,  physical 
mechanisms  linking  sunspots  and  terrestrial  cyclones  have 
eluded  investigators  and  although  Kullmer  (1933,  1943) 
updated  his  analyses  of  cyclone  frequencies,  we  are  no 
nearer  to  the  predictive  capacity  he  envisioned. 

Interest  in  and  appreciation  of  the  significance 
of  the  climatology  of  cyclones  was  revived  during  and 
after  World  War  II  (Petterssen,  1941?  Weightman,  1945; 
Miller,  1946;  Miller  and  Mantis,  1947;  Saucier,  1949; 

James,  1952;  Hurley,  1954) .  This  work  culminated  with 
Klein's  (1957,  1958)  detailed  studies  of  spatial  distri¬ 
bution  of  cyclone  frequency  and  Hosier  and  Gamage ' s  (1956) 
determination  of  the  secular  history  (1905-1954)  of  cyclone  frequencies. 

During  the  1960s  attention  shifted  to  the  marine 
sector.  Mather,  Adams  and  Yoshioka  (1964)  and  Mather, 

Field  and  Yoshioka  (1967)  summarized  the  climatology 

of  damaging  Atlantic  coast  cyclones  and  reported  a  secular 

change  toward  more  frequent  damaging  storms 

caused  by  lowering  of  the  average  central  pressure  of 

storms  off  the  Atlantic  coast.  More  recently,  Zishka 

and  Smith  (1980)  reported  a  general  trend  of  lower  central 
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pressure  for  cyclones  over  the  North  American  sector. 

Hayden  (1975) ,  Resio  and  Hayden  (1975) ,  and  Hayden  (1980a) 
confirmed  the  trend  over  the  last  five  decades  toward 
more  frequent  and  more  intense  cyclones  off  the  U.S.  Atlantic 
coast.  For  the  Atlantic  coastal  region  a  shift  in  the 
mean  path  of  movement  of  cyclones  has  been  documented 
(Resio  and  Hayden,  1975;  Dickson  and  Namias,  1976;  Hayden, 
1980a)  . 

In  the  last  few  years  improved  climatological 
assessments  of  cyclone  frequencies  have  entailed  more 
extensive  and  intensive  grids  for  sampling  (Reitan,  1974, 
1979;  Zishka  and  Smith,  1980)  and  more  extensive  historical 
coverage  of  cyclone  frequency  data  fields  (Hayden, 
1980a) . 

b.  Climate  Prediction 

The  predictability  of  climate  is  conditional  on 
establishing  and  having  foreknowledge  of  forcing  pro¬ 
cesses,  e.g.  sunspots  (Kullmer,  1943) ,  chandler  motions 
of  the  Earth's  axis  of  rotation  (Bryson  anr*  Starr,  1977),  sea 
surface  temperature  fields  and  atmospheric  circulation 
(Namias,  1974),  analog  climate  states  (Barnett  and  Preisen- 
dorfer,  1978)  or  establishing  the  existence  of  seasonal  or 
longer  climatological  persistence  such  as  identified  by  Sir 
Gilbert  Walker  and  E.W.  Bliss  (1932)  and  erployed  in  the  present  study. 
Resolving  the  causes  of  climatological  persistence  and  its 
occasional  failures  returns  one  to  the  issue  of  climate¬ 
forcing  processes. 


c.  Multivariate  Statistical  Prediction  Models 


Early  use  of  multivariate  techniques  in  prediction 
have  focused  on  prediction  of  prehistoric  atmospheric 
conditions  from  modern  data.  Fritts  et  al.  (1971)  used 
canonical  correlations  between  patterns  of  hemispheric 
pressure  and  patterns  of  tree-ring  growth  to  extend 
(synthesize)  backward  pentad  mean  seasonal  pressure  charts 
to  1700.  In  similar  fashion  Webb  and  Bryson  (1972)  used 
multivariate  statistics  to  reconstruct  precipitation, 
temperature  and  snowfall  conditions  during  the  Holocene 
in  the  upper  Midwest  from  modern  pollen  and  weather  data 
and  from  fossil  pollen  stratigraphies.  Similar  procedures 
have  been  employed  to  predict  (hindcast)  sea  surface 
temperatures  at  the  end  of  the  last  ice  age  (18,000  BP) 
from  plankton  in  sea  floor  cores  (CLIMAP,  1976) . 

Torranin  (1972)  used  canonical  correlation  to  predict 
hydrological  conditions  in  the  mountains  of  California 
from  vector  patterns  of  sea  surface  temperature  in  the  North 
Pacific.  Gilman  (1957)  used  principal  component  analyses 
of  pressure,  precipitation  and  temperature  fields  to  predict 
one  variable  from  the  others.  In  the  current  study, 
correlations  between  principal  components  of  winter  and 
summer  cyclone  frequency  patterns  define  three  modes  of 
climatological  persistence  such  that  knowledge  of  the 
vector  patterns  of  cyclone  frequencies  in  one  season  permits 
estimation  of  frequency  patterns  in  the  following  season. 
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3 .  The  Data 

Cyclone  frequencies  for  winter  (October  -  March)  and 
summer  (April  -  September)  for  the  years  1885-1979  were 
tabulated  for  the  74,  2 . 5°-latitude  by  5 . 0°-longitude 
grid  cells  comprising  the  study  areas  shown  in  Figure  1. 
From  monthly  charts  of  the  "Tracks  of  the  Centers  of 
Cyclones  at  Sea  Level"  published  by  the  Monthly  Weather 
Review  and  in  recent  years  by  the  Mariners  Weather  Log, 
seasonal  totals  of  cyclones  passing  through  each  grid 
cell  were  recorded.  Multiple  entries  of  a  given  storm  in 
a  grid  cell  were  ignored.  The  resulting  data  matrices 
have  the  dimensions  of  74  grid-cell  variables  and  94 
yearly  cases.  Frequencies  were  not  adjusted  for  latitude 
variations  in  grid-cell  area  because  of  distortions 
involved  in  such  adjustments  (Hayden,  1980b) . 

4 .  Long-Term  Means  and  Variances 

Figure  2  shows  the  mean  winter  and  summer  frequency 
of  cyclones  per  grid-cell  for  the  period  1885-1979.  The 
winter  and  summer  patterns  are  similar  in  all  respects 
except  in  magnitude.  Cyclones  are  more  frequent  in  winter. 
These  seasonal  mean  frequency  patterns  are  similar  to  the 
annual  means  shown  by  Hayden  (1980a) .  In  general  cyclones 
increase  with  latitude  and  proximity  to  the  east  coast 
where  a  dominant  frequency  maximum  axis  is  found.  The 
location  of  this  frequency  maximum  is  perhaps  associated 
with  the  dynamic  affects  of  the  Appalachian  Mountains  and 
baroclinic  affects  of  the  coastal  zone. 
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The  standard  deviation  of  winter  and  summer  cyclone 
frequencies  is  presented  in  Figure  3.  Like  the  patterns 
of  the  means,  the  standard  deviations  for  winter  and  summer 
are  similar  in  form  and  bear  great  similarity  to  the  annual 
pattern  published  earlier  (Hayden,  1980a)  .  The  dominant 
feature  of  the  chart  is  the  pronounced  frequency 
maximum  off  the  Atlantic  coast.  The  axis  of  this  maximum 
is  seaward  of  the  maximum  in  the  mean  field  (Fig.  2)  and 
apparently  rotated  clockwise  several  degrees.  The 
location  of  the  standard  deviation  maxima 
closely  approximates  the  locus  of  the  Gulf  Stream. 

The  geographic  relationship  between  the  Gul-d  Stream  and  the 
variance  maxima  in  cyclone  frequencies  is  interesting  in 
view  of  the  season-to-season  persistences  detailed  in  sub¬ 
sequent  sections.  The  standard  deviations  landward  of  the 
east  coast  are  rather  modest  suggesting  that  the  axis  of 
storm  occurrences  may  expand  eastward,  but  in  general  not 
westward,  of  the  Appalachians. 

5.  Principal  Component  Analyses  of  Winter  and  Summer 
Frequency  Data 
a.  The  Analysis 

Principal  component  analyses  have  successfully 
resolved  the  variance  structure  in  multivariate,  geo¬ 
physical  data  (Kutzbach,  1967;  Fritts  et  al.,  1971;  Resio 
and  Hayden,  1975;  Hayden,  1980a) ,  and  in  terms  of  least 
square  errors,  this  type  of  analysis  provides  a  method  for 


determining  patterns  in  large  data  fields  (Lorenz,  1956; 
Gilman,  1957;  Kutzbach,  1967) .  The  objective  of  the  analysis 
is  to  isolate  characteristic,  recurrent,  and  independent 
modes  of  covariance  among  variables  into  a  new  set  of 
independent  variables.  Basically,  the  analysis  transforms 
a  set  of  intercorrelated  variables  into  a  new  coordinate 
system  in  which  the  axes  are  linear  combinations  of  the 
original  variates  and  are  mutually  orthogonal. 

To  prevent  those  grid  cells  with  high-mean  cyclone  freq¬ 
uencies  (high  latitudes)  from  dominating  the  total  variance 
and  consequently  from  dominating  the  eigenvector  forms,  the 
correlation  matrix  was  used  rather  than  the  covariance  matrix. 
The  procedures  for  calculation  of  the  eigenvectors  follow 
those  of  Kutzbach  (1967)  and  Vincent  et  al.  (1976). 

Principal  component  analysis  provides  a  description 
of  the  uncorrelated  major  modes  of  variation  in  the  data  set.  Typically, 
each  component  is  identified  with  some  property  of  the  data 
field.  The  analysis  also  provides  an  index  which  measures  the 
importance  of  each  component  within  each  case.  Finally  the 
analysis  provides  an  estimate  of  the  total  percent  of  variance 
in  the  data  set  which  can  be  explained  on  the  basis  of  each 
component.  Separate  principal  component  analyses  were  run  for  the 
winter  and  summer  data  matrices, 
b.  Significance  of  the  Eigenvectors 

Following  the  lead  of  Barnett  and  Preisendorfer  (1978) , 

94  year-by-74  grid  cell  matrices  of  randomly  selected  (winter 
and  summer)  cyclone  frequencies  were  prepared  and  subjected  to 
eigenvector  analyses.  The  percentage  of  variance  explained 
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by  each  of  the  first  20  eigenvectors  calculated  as 
compared  with  the  percentages  of  variance  explained 
by  the  eigenvectors  of  actual  winter  and  summer  cyclone 
frequency  occurrences  is  shewn  in  Figure  4.  For  both  winter 
.J  and  summer  analyses  the  first  5  eigenvectors  accounted 

!  for  more  variance  than  the  first  six  "random"  cyclone 

^  eigenvectors.  Only  the  first  four  eigenvectors  are 

i 

■*i  studied  in  this  report. 

^  c.  The  Eigenvectors 

The  percentage  of  variance  and  the  cumulative 
percentage  of  variance  explained  by  the  first  four 

it 

!  eigenvectors  for  the  winter  and  summer  seasons  is 

4  given  in  Table  1.  As  indicated,  the  first  four 

winter  and  summer  eigenvectors  account  for  52  percent 
and  49  percent  of  the  total  variance,  respectively. 

The  problem  has  thus  been  reduced  from  two  74-variable 
problems  in  which  the  variables  are  intercorrelated, 
to  two  4-variable  problems  in  which  each  new  variable 
is  orthogonal  and  thus  statistically  independent. 
Higher  order  eigenvectors  explained  small  fractions 
of  the  total  variance  and  are  not  used  in  the  model 
developed  here. 

The  winter  and  summer  eigenvectors  corres¬ 
ponding  to  the  eigenvalues  .1  through  4  are  mapped 
in  Figures  5,  6,  7,  and  8.  '’’he  first  eigenvectors 

of  winter  and  summer  analyses  (Fig.  ?)  are 


similar  in  form  and  amplitude  and  indicate  that 
the  dominant  mode  of  cyclone  frequency  pattern 
variation  about  the  mean  is,  with  positive  weight¬ 
ings,  an  abundance  of  cyclones  over  marine  areas 
occurring  at  the  expense  of  cyclone  frequencies 
over  the  continent.  Years  with  negative  weight¬ 
ings,  in  both  winter  and  summer,  would  be  character 
ized  by  fewer  than  average  cyclones  over  marine 
areas  and  above  average  frequencies  over  the 
continent.  Using  annual  data  it  was  shown  that 
this  pattern  is  not  the  result  of  poor  cyclone 
identification  over  marine  areas  (Hayden, 

1980a) . 

The  second  eigenvectors  of  the  winter  and 
summer  data  matrices  (Fig.  6)  have  positive  values 
over  the  entire  field  and  as  such  indicate,  for 
positively-weighted  years,  a  general  increase  in 
cyclone  numbers  and  in  negative ly-weighted  years, 
fewer  than  average  cyclones.  The  dominant  feature 
the  maps  of  the  second  eigenvectors  is  the  maximum 
centered  over  the  coast  of  the  mid  -  Atlantic. 
This  pattern:  is  rather  similar  to 
the  chart  of  cyclogenesis  frequency 
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published  by  Petterssen  (1941) .  As  was  the  case  with  the 
first  pair  of  eigenvectors,  the  second  eigenvectors  are 
essentially  the  same  in  form  and  amplitude.  There  is  no 
indication  of  seasonal  differentiation. 

The  third  eigenvector  in  both  winter  and  summer  (Fig.  7) 
contrasts  storms  tracking  eastward  from  the  Great  Plains 
up  the  Ohio  Valley  and  exiting  the  coast  between  Cape 
Hatteras  and  Cape  Cod,  with  frequent  cyclones  in  the  Gulf 
coast  region.  Again  no  seasonal  differentiation  is  evident. 

The  fourth  eigenvectors  (Fig.  8)  contrast  the  frequent 
cyclone  occurrence  in  the  Great  Lakes-Saint  Lawrence  River 
region  with  cyclones  out  of  the  Great  Plains.  In  years 
with  positive  weightings,  storms  are  more  frequent  than 
average  in  southern  Canada.  Negatively-weighted  years  show 
an  increase  in  storms  along  the  Colorado  storm  track. 

The  first  four  eigenvectors  of  the  winter  and  summer 
cyclone  frequency  matrices  constitute  four  pairs  of  new, 
orthogonal  variables  that  account  for  about  50  percent  of 
the  variance  in  the  original  data  sets.  These  new  variable 
pairs  are  functionally  similar  in  form  and  amplitude.  There 
is  no  evidence  of  seasonal  differences  in  the  patterns. 

The  four  vector  pairs  closely  resemble  vectors  calculated 
from  an  annual  data  matrix  published  earlier  (Hayden,  1980a) . 
d.  The  Time  Series 

The  time  series  weightings  of  each  of  the  four  eigen¬ 
vector  pairs  for  each  year  (1885-1979)  are  shown  in  Figure  9. 
The  annual  weightings  of  the  first  eigenvectors  in  winter 
and  summer  show  broad,  quasi-sinusoidal  variation  of  period 
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slightly  greater  than  100  years.  This  variation  suggests 
that  in  both  winter  and  summer,  cyclone  frequencies 
declined  over  marine  areas  and  increased  over  the  continent 
between  1885  and  1925.  Since  1925,  cyclone  frequencies 
over  marine  areas  have  increased  and  frequencies  over  the 
continent  have  declined.  Visual  inspection  of  Fig.  8A 
for  winter  and  summer  indicates  covariations  at  both  high 
and  low  frequencies. 

Annual  weightings  of  the  second  winter  and  summer 
eigenvectors  (Fig.  9  )  also  exhibit  century-long  scales  of 
variation  from  negative  values  at  the  beginning  of  the 
time  series  to  positive  values  after  1910  and  back  to 
negative  values  around  1960.  If  the  interpretation  of  the 
second  eigenvectors  as  representing  coastal  cyclogenesis 
is  correct,  then  there  should  have  been  an  increase  in 
coastal  cyclogenesis  during  the  first  five  decades  of 
the  20th  century.  This  observation  is  supported  by  Resio 
and  Hayden  (1975)  who  show  an  increase  in  cyclogenesis  along 
the  mid- Atlantic  coast.  The  second,  eigenvector  time  series  for  both  winter 
and  surrmer  ("Ig.  9)  also  shows  remarkable  similarities  to  Reitan's  (1974) 
tabulations  of  total  cyclone  numbers  for  the  western  part  of  the  Northern 
Hemisphere.  As  was  the  case  with  the  first  eigenvector  pairs,  both  high 

and  low  frequenc  '  covariations  are  evident. 

The  time  series  weightings  for  third  and  fourth 
eigenvector  pairs  (Fig.  9)  also  show  winter-summer 
similarities.  To  quantify  the  similarities  between  the  first 
four  eigenvector  weightings  for  the  winter  and  summer  time  series, 
correlation  coefficients  were  calculated  for  both  winters 


versus  the  following  summers  and  for  summers  versus  the 
following  winters.  These  correlation  coefficients  are 
summarized  in  Table  2.  Clearly  there  are  at  least  three 
synoptic  scale  modes  of  seasonal  persistence  evident  from 
these  analyses.  Given  the  essential  identity  between  the 
winter  and  summer  eigenvector  pairs  and  the  significant 
correlation  between  their  respective  time  series,  a  pre¬ 
dictive  model  can  be  constructed  and  evaluated  using  the 
first,  second  and  fourth  eigenvectors. 


6.  The  University  of  Virginia  Climate  Prediction  Model 
a.  The  Predictive  Equations 

The  cyclone  frequency  values  at  each  of  the  74  grid  cells 
(Fig.  1)  for  a  given  winter  (Cw)  may  be  estimated  from 
the  following: 


_  *  *  * 

C  =  X  +  a.,  a  E,  +  a_  a  E_  +  a.  ,  cr  E. 

w  w  lw  w  lw  2w  w  2w  4w  w  4w 


where  X  is  the  matrix  of  74  long-term  (1885-1979)  mean 
w 

winter  cyclone  frequencies  (Fig.  2) ;  is  the  standard 
deviation  matrix  of  X  (Fig.  3);  E.  are  the  four  winter 

1W 

★ 

eigenvectors  (Figs.  5,  6,  7,  and  8);  and  a^w  are  the  predicted 

* 

eigenvector  weightings.  The  a^w  are  calculated  from 
regressions  between  the  winter  and  summer  eigenvector 
weighting  time  series  (Fig.  9)  and  the  eigenvector  weightings  for 
the  sumner  preceding  the  winter  forecasted.  A  similar  equation  is 


structured  for  the  prediction  of  summer  cyclone  frequency 
patterns.  In  short  this  forecast  method  states  that  the 
cyclone  frequency  pattern  for  the  coming  season  should  be 
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similiar  to  the  average  for  that  season  with  the  addition 
of  three  persistence  components  based  upon  the  occurrence 
of  cyclones  during  the  preceding  season.  With  the  Univer¬ 
sity  of  Virginia  Climate  Prediction  Model,  a  forecast  of 
cyclone  frequencies  is  produced  for  each  2 . 5 “-latitude 
by  5 . 0“-longitude  grid  cell  of  the  study  area  (Fig.  1). 

The  cyclone  frequency  forecast  for  each  grid  cell  can  be 
compared  to  the  actual  cyclone  occurrence  and  forecast 
skill  can  be  quantified, 
b.  Model  Evaluation 

In  order  to  evaluate  the  predictive  utility  of  the 
University  of  Virginia  Climate  Prediction  Model,  test 
forecasts  were  made  using  independent  data.  From  the 
total  of  94  years  of  winter  and  summer  data  available, 
every  ten  years  were  removed  from  the  data  matrix.  The 
new  smaller  winter  and  summer  matrices  of  74  variables  and 
84  cases  were  subject  to  principal  component  analyses. 

The  eigenvectors  derived  and  the  weightings  for  each  season 
were  calculated  and  regressed  against  the  previous  season. 
The  eigenvector  weightings  for  summer  were  used  to  estimate 
the  weightings  for  the  following  winter  and  vice  versa. 

The  forecasted  weightings  were  then  applied  to  the  pre¬ 
dictive  equations  and  a  cyclone  frequency  forecast  was 
produced.  The  ten  independent  years  of  data  were  then 
returned  to  the  matrix  and  a  new  set  of  independent  years 
was  removed,  the  model  was  rebuilt,  and  a  second  forecast 
was  produced.  This  was  repeated  until  all  94  years  of 


record  were  predicted  for  both  winter  and  summer.  While 
this  method  of  selection  of  independent  data  is  an  adequate 
test  it  differs  from  prediction  of  the  next  ten  years  (1980s)  . 
To  approximate  this  forecast  situation,  the  1970s  were 
removed  from  both  the  winter  and  summer  data  matrices . 

The  models  and  predictive  equations  were  reformed  and  each 
year  was  forecast  in  turn.  On  completion  of  these  forecast 
exercises,  the  results  were  compared  to  historical  occurrences 
and  forecast  skill  was  calculated, 
c.  Predictive  Skill 

Numerous  methods  have  been  advanced  to  quantify  estimates 
of  forecast  skill  (Brier  and  Allen,  1951;  Vernon,  1953) 
and  as  noted  by  Brier  and  Allen  the  method  selected  depends 
on  the  purpose  of  verification.  The  purpose  here  is  twofold. 
First,  it  is  essential  to  establish  the  relative  level  of 
reliability  of  the  University  of  Virginia  Climate  Forecast 
Model  and  secondly,  identify  forecast  errors  and  determine 
their  nature  and  possible  cause.  The  basic  test  of  forecast 
skill  employed  here  is  the  percent  of  correct  forecasts 
(percent  skill  score)  relative  to  the  use  of  the  climatological 
mean  as  a  forecast.  Thus  while  the  model  provides  a  numerical 
value  of  cyclone  frequency,  the  skill  test  used  here  evaluates 
only  the  accuracy  of  forecasts  of  the  form  greater  or  less 
than  the  mean.  For  comparison,  Heidke  skill  scores  , 
absolute  errors  of  the  forecast,  root  mean  square  errors  of 
the  forecast,  deviation  skill  score  and  quadratic  skill 
scores  were  also  calculated  and  will  be  summarized. 


For  each  season  of  each  of  the  94  years  forecasted 
the  forecast  cyclone  frequency  at  each  grid  cell  was 
compared  to  the  actual  occurrence  of  cyclones  and  the  long¬ 
term  mean  frequency.  The  percentage  of  the  74  grid  cells 
was  calculated  for  each  season  and  year  where  the  model 
provides  a  better  estimate  of  actual  cyclone  occurrences 
than  the  mean.  This  aggregate  skill  is  referred  to  here 
as  global  skill  following  Barnett  and  Preisendorfer  (1978) . 
The  average  (1885-1979)  percent  skill  score  was  also 
determined  for  each  grid  cell  for  each  season  and  is 
referred  to  as  local  skill, 
d.  Global  Skill 

Figure  10  gives  the  grid  average  (global) ,  percent  skill 
scores  for  season  by  each  year.  The  average  percent  skill 
score  for  the  94  years  is  66  percent  for  both  winter  and 
summer  season  predictions.  Frequency  histograms  of  yearly 
global  skill  scores  by  5  percent  class  intervals  are  given 
in  Figure  11.  The  distribution  of  percent  skills  is  near 
normal.  For  winter  forecasts,  92  percent  of  the  years 
forecasted  had  positive  skill  while  87  percent  of  the  summer 
forecasts  had  positive  skill.  All  years  with  negative 
global  forecast  skill  were  isolated  for  additional  study 
(skill  was  calculated  for  only  the  first  half  of  the  season 
forecasted) .  Half  of  the  failure  years  were  found  to  have 
positive  forecast  skill  for  the  first  three  months  and 
percentage  improvements  in  forecast  skill  were  achieved  for 
92  percent  of  the  original  forecast  failures. 
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The  population  of  forecast  failure  years  was  also 
subjected  to  a  new  model  construction.  In  this 
exercise  only  the  actual  cyclone  occurrences  in  the 
latter  half  of  the  antecedent  season  were  used  to 
derive  the  forecast  of  the  subsequent  season.  Again 
about  half  of  the  failure  years  re  -  forecasted 
showed  positive  skill.  These  tests  indicated  that 
major,  mid-season  change  in  the  general  circulation 
cause  forecast  failure  in  the  subsequent  season  and 
that  general  circulation  changes  within  a  season  are 
a  frequent  cause  of  failure  to  achieve  positive  forecast 
skill  for  that  season. 

The  average  percent  skill  for  winter  and  summer 
as  compared  with  Heidke  skill  scores  (Heidke,  1926) , 
absolute  errors  and  root  mean  square  errors  (RMSE) , 
deviation  skill  and  quadratic  score  (Vernon,  1951) 
are  given  in  Table  3.  As  noted  by  Serf ling  (1975) 
"statistical  models  for  meteorological  phenomena 
typically  perform  poorly  in  extremes."  When  all  94 
years  forecasted  are  considered, that  is  the  case 
here.  The  quadratic  and  deviation  skill  scores 
penalize  failure  to  predict  extremes.  With  the  deviation 
skill  score,  the  penality  of  incorrect  forecasts  increases 
linearly  with  increasing  departure  from  the  mean  and  for 
the  quadratic  skill  score  the  penality  increases  as  the 
square  of  the  departure  from  the  mean.  Figure  12 
illustrates  the  relationship  between  global  percent 
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skill  scores  and  quadratic  skill  scores  for  the  94 
stammer  predictions  and  the  9  4  winter  predictions. 

When  both  the  percent  and  quadratic  skill  scores  are 
positive,  the  two  measures  are  linearly  related. 

There  is  no  apparent  relationship  between  percent 
and  quadratic  skill  score  when  only  cases  with 
negative  quadratic  skill  scores  are  considered.  It 
seems  clear  that  the  model  being  tested  here  does  a 
good  job,  in  general,  predicting  the  form  or  pattern 
of  departure  from  the  mean  and  when  forecast  skill  is 
positive  a  good  job  of  predicting  magnitude  of  the 
departures  from  long-term  averages. 

While  the  procedure  of  selecting  test  years  from 
within  the  data  record  is  an  independent  test,  questions 
still  remain  about  extension  of  the  forecast  into  the 
future.  With  the  1970s  removed  from  the  data  record 
the  model  was  rebuilt  and  the  winter  of  1969-1970 
and  summer  of  1970  were  predicted  and  percent  skill 
scores  were  calculated.  The  1970  year  was  then  added 
back  to  the  data  matrix,  the  model  was  rebuilt  and 
the  next  winter  and  summer  were  predicted.  Each  year 
of  the  1970s  was  computed  in  this  manner  in  an  effort 
to  mimic  a  real  forecast  situation.  In  Figure  13  the 
resulting  global  percent  skill  scores  for  each  year 
are  compared  with  those  given  in  Figure  10.  Heidke 
deviation  and  quadratic  skill  scores  are  also  shown. 

In  both  winter  and  summer  the  forecast  skill  is  about 
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the  same  as  that  calculated  from  within  the  data  record. 

The  model  thus  is  equally  good  in  forecasting  into  the 
future.  Figure  14  compares  the  global  percent  skill 
scores  for  the  forecast  of  the  1970s  with  the  quadratic, 
deviation  and  Heidke  global  skill  scores.  The  various 
scores  parallel  each  other.  This  seems  consistent  with 
the  data  shown  in  Figure  12.  Figure  12  clearly  suggests 
that  when  the  global  percent  skill  score  is  high  both 
the  global  quadratic  and  deviation  skill  scores  are  also 
high.  When  global  percent  skill  score  falls  below  50  percent 
the  quadratic  and  deviation  skill  scores  become  negative. 
These  results  suggest  that,  contrary  to  Serf ling's  (1975) 
dictum  about  statistical  models'  general  failure  to  predict 
extremes,  the  University  of  Virginia  Climate  Prediction 
Model  does  a  good  job  of  predicting  extreme  values  when 
a  positive  percent  skill  is  achieved. 

e.  Local  Skill 

Figure  15  shows  the  geographic  variation  (local  skill) 
of  forecast  skill  for  winter  and  summer.  In  both  charts, 
positive  skill  is  found  for  all  grid  cell  locations.  (Histo¬ 
grams  of  winter  and  summer  grid-cell  percent  skill  scores 
by  5  percent  class  intervals  are  given  in  Figure  15.)  In 
general,  winter  forecasts  have  a  greater  range  of  percent 
skill  about  the  chart  mean  of  66  percent  than  does  summer. 

The  isopleth  field  for  winter  is  also  more  highly  organized 
than  the  summer  field, 
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f.  Long-term  versus  short-term  global  skill 

The  time  series  of  eigenvector  weightings  (Fig.  9) 
exhibit  long-term  trends,  i.e.  year-to-year  and  longer 
persistences.  Inspection  of  Figure  9  also  reveals 
numerous  winter- summer  synchronous  departures  from  the 
general  trends.  Third-order  polynomial  curves  were 
fitted  to  the  tine  series  of  eigenvector  weightings  . 
using  a  least  square  criterion.  The  polynomial  curves 
were  then  subtracted  from  the  original  time  series. 

The  time  series  of  polynomial  departures  (PD)  were 
correlated  with  the  time  series  of  global  skill  (S)  for 
the  following  season.  These  correlations  are  summarized 
in  Table  4.  with  negative  correlations,  global  skill 
declines  with  increasing  departures  from  the  polynomial 
trend  and  indicates  a  dominance  of  long-term  persistence 
in  the  predictions.  With  positive  correlations,  global 
skill  increases  with  increasing  departures  from  the  trends 
indicating  high  frequency  (short-term)  contributions  to 
overall  skill.  Table  4  indicates,  however,  that  strong 
positive  and  negative  correlations  are  not  found  and  thus 
there  is  no  dominance  of  high  or  low  frequency  persistences 
in  the  skill  achieved.  Fven  if  the  correlations  were 
more  robust  and  significant,  the  slopes  on  the  regression 
are  modest  (Table  4).  For  example,  use  of  the  polynomial 
curve  of  eigenvector  number  1  weightings  rather  than  the 
actual  weightings  would  net  a  potential  skill  improvement 
of  about  2.6  percentage  points.  In  similar  fashion,  using 
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just  the  departures  from  the  polynomial  trend  for  eigenvector 
4  would  net  a  potential  skill  improvement  of  4  percentage  points. 
Unfortunately  the  correlations  are  weak  and  we  must  conclude  that 
the  skill  achieved  with  the  University  of  Virginia  Climate 
Prediction  Model  is  the  result  of  both  long-  and  short-term 
persistences  in  the  eigenvector  weightings. 

g.  Forecast  Skill  Relative  to  Chance,  Persistence  and  the  Previous 
Season  as  a  forecast. 

In  earlier  sections  skill  was  assessed  relative  to  the  long¬ 
term  mean  (1885-1979) .  Here,  the  skill  achieved  relative  to 
chance,  persistence,  and  the  previous  season  as  forecasts  is 
evaluated.  In  each  of  these  computations  of  forecast  skill  a 
94  year-by-74  grid  cell  matrix  of  winter  and  summer  cyclone 
frequency  was  prepared  as  a  base  of  comparison.  The  cyclone 
frequency  matrix  of  random  (chance)  values  was  discussed  in 
an  earlier  section. 

The  matrices  for  the  test  of  skill  relative  to  persistence 
were  constructed  from  the  actual  data.  Third-order  polynomials 
were  fitted  to  the  94  year  time  series  for  each  of  the  74  grid 
cells  for  both  the  winter  and  summer  data.  From  the  148  poly¬ 
nomial  curves,  values  for  each  grid  cell  for  each  year  and  season 
were  extracted.  The  resulting  winter  and  summer  matrices  were 
used  in  the  persistence  tests  for  model  skill. 

The  simplest  forecast  to  test  model  skill  against  is  that 
the  coming  season  will  be  like  the  previous  season.  Unfortunately 
cyclone  frequencies  are  higher  in  all  grid  cells  in  winter  than 
summer.  To  make  adjustments  for  this  seasonal  difference,  the 
winter  matrix  was  multiplied  by  the  ratio  of  the  summer  to  winter 
long-term  means  for  each  grid  cell.  The  summer  matrix  was 
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multiplied  by  the  ratio  of  winter  to  summer  long-term  means. 

In  short,  the  winter  of  1978-79,  when  adjusted  to  summer  frequency 
levels,  becomes  the  forecast  for  the  summer  of  1979. 

The  University  of  Virginia  Climate  Prediction  Model  achieves 
a  higher  level  of  skill  than  the  mean,  chance,  persistence  or 
the  previous  season  as  a  forecast  (Table  5) .  Skill  relative  to 
persistence  shows  the  most  modest  gains.  This  is  not  surprising 
because  the  model  is  basically  a  persistence  forecast  system. 

Based  on  the  results  summarized  in  Table  5,  it  apppears  that  about 
two-thirds  of  the  skill  is  due  to  long-term  persistence  and  about 
one-third  to  shorter-term  persistence.  The  previous  season, 
however,  is  not  a  very  good  forecast  of  the  season  which  follows. 
Lag  regressions  of  the  time  series  shown  in  Figure  9  suggest 
that  model  skill  for  the  time  interval  near  or  slightly  longer 
than  1  year  may  be  better  than  the  6-month  seasons  studied  here. 
The  deviation  and  quadratic  skill  scores  are  negative  for  the 
model  relative  to  persistence.  This  results  because  of  the  higher 
number  of  failure  years  using  this  most  stringent  base  of  compari¬ 
sons  (see  Figure  12) .  Most  of  the  years  with  a  percent  skill 
score  of  60  or  higher  for  the  test  of  skill  relative  to  per¬ 
sistence  had  positive  deviation  and  quadratic  skill  scores. 

The  only  winter/summer  disparity  in  skill  scores  comes 
from  comparison  of  model  skill  achieved  with  the  previous 
season  as  a  prediction.  Apparently  summers  are  not  as  good 
as  estimators  of  the  following  winter  as  visa  versa.  Or  put 
more  directlv,  the  character  of  winter  more  frequently  extends 
over  to  summer  than  does  summer  to  winter. 


h.  The  Winter  of  1979-1990 

The  first  "future"  forecast  made  using  the  Uni¬ 
versity  of  Virginia  Climate  Prediction  Model  was  for 
the  winter  of  1979-1990.  The  forecast  was  completed 
in  November  of  1979  and  has  been  verified.  Percent, 
Heidke,  deviation  and  quadratic  skill  scores  were  94  per 
cent  .68,  .18,  and  .33  respectively.  The  reduction 
in  the  root  mean  square  error  of  the  forecast,  relative 
to  the  mean  as  a  forecast,  was  .56  cyclones  per 
grid  cell.  The  reduction  in  the  absolute  error  of  the 
forecast  relative  to  the  mean  was  0.44  cyclones  per 
grid  cell.  Figure  17  shows  comparisons  of  the  1979- 
19  80  predicted  and  actuaJ.  cyclone  frequencies,  the 
difference  between  the  predicted  and  the  long-term  mean 
and  the  difference  between  the  actual  and  mean  fre¬ 
quencies,  and  the  difference  between  the  actual  occur¬ 
rences  and  the  predicted  occurrences.  Forecast  skills 
are  also  listed  in  Figure  17.  Twelve  of  the  74  grid 
cells  were  not  forecasted  correctly  (Figure  14F) . 

For  example,  at  37.5°N,  97.58W  fewer  storms  than  normal 
were  predicted  and  fewer  than  normal  occurred;  the 
cyclone  frequencies  at  this  location  were  4.74  fewer 
than  predicted.  The  anomaly  was  under-estimated 
at  37.5°N,  67.5°W  where  more  than  average  cyclone 
numbers  were  predicted,  and  more  than  normal 
occurred.  The  anomaly  at  this  site  was  also 
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underestimated .  Overall  forecast  skill  for  the  1979-1980 
winter  was  very  good. 

7.  Discussion  and  Conclusion 

The  University  of  Virginia  Climate  Prediction  Model 
provides  a  forecast  of  cyclone  frequency  patterns  for 
eastern  North  America  and  the  western  North  Atlantic.  For 
both  winter  and  summer,  area  average  percent  forecast  skill 
scores  are  66  percent.  Geographically  the  model  shows 
positive  forecast  skill  in  all  areas  and  as  expected 
forecast  skill  varies  in  both  space  and  time.  The  model 
essentially  produces  a  prognostic  chart  of  the  expected 
spatial  cyclone  frequency  distribution  (Fig.  17) .  The 
interpretation  of  these  charts  in  terms  of  warmer  and  cooler 
and  wetter  and  drier  regions  requires  additional  study. 
Schematic  models  such  as  those  suggested  by  Klein  (1948) 
should  prove  particularly  useful  as  chart  interpretation 
aids.  In  addition,  statistical  pattern  analysis  methods 
(Biasing,  1975;  Kutzbach,  1967;  Fritts  et  al. ,  1971)  and  robust 
data  collections  should  permit  the  product  of  the  University 
of  Virginia  Climate  Prediction  Model  to  derive  second  order 
models.  For  example,  preliminary  study  indicates  that  the 
yearly  weightings  of  the  first  eigenvector  of  cyclone 
frequencies  are  highly  correlated  with  heating  degree  day 
departures  from  the  historical  trend  line  for  Washington, 
D.C.,  and  is  also  highly  correlated  with  the  frequency  of 
one  inch  or  greater  snowfalls  in  Richmond,  Virginia. 
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Considerable  study  is  required  to  evaluate  the  reliability 
and  value  of  such  derived  forecasts. 

The  University  of  Virginia  Climate  Prediction  Model 
is  a  persistence  forecast  system.  During  the  early  phases 
of  this  investigation  the  discovery  of  season-to-season 
components  of  persistence  in  the  fields  of  cyclone  frequencies 
was  rather  surprising.  Aside  from  the  early  work  of  Sir  Gilbert 
Walker  and  E.  W.  Bliss  (1932)  anda  few  scattered  reports,  the  liter¬ 
ature  is  depauperate  in  citations  on  seasonal  persistences 
of  atmospheric  properties.  Those  seasonal  persistences 
which  are  reported  are  generally  of  synoptic  scale  or 
larger.  Namias  (1960)  refers  to  "contrasting  'regimes' 
on  the  order  of  a  month,  in  which  weather  processes 
persistently  recur,"  and  Ballenzweig  (1959)  notes  that 
over  a  period  of  a  month  or  a  season,  "it  is  frequently 
found  that  the  tracks  of  tropical  cyclones  cluster  about 
certain  preferred  axes."  Namias,  Dunn  and  Simpson  (1955) 
refer  to  this  cluster  about  a  preferred  axis  as  a  "closely 
knit  family."  It  is  clear  that  persistence  in  the  move¬ 
ment  of  cyclones  within  a  season  has  been  observed  and  this 
is  visually  born  out  in  inspections  of  the  monthly  charts 
of  the  centers  of  cyclones.  What  is  less  clear  and  more 
important  is  the  physical  basis  for  the  observed  persistence. 

Some  clues  as  to  the  physical  basis  for  cyclone 


frequency  pattern  persistence  can  be  derived  from  the  charts 
of  means  and  variances  (Figs.  2  and  3).  In  Figure  2,  a 
pronounced  frequency  minima  is  evident  along  the  Appalachian 


Mountains.  Colucci  (1976)  in  a  high  spatial  resolution 
cyclone  study  of  the  region  confirms  this  frequency  minima. 
Colucci  also  notes,  however,  that  cyclone  deepening  over 
the  Appalachians  is  also  typical.  It  would  appear  then 
that  the  Appalachians  exert  dynamic,  albeit  temporally 
constant,  effects  on  cyclone  frequency  patterns.  The 
maximum  frequencies  in  the  mean  field  (Fig.  2)  closely 
parallel  the  coast,  also  a  static  feature.  The  maximum 
in  the  variance  field  (Fig.  3) ,  however,  is  not  congruent 
with  the  means.  The  variance  maximum  closely  parallels 
the  Gulf  Stream.  Colucci  (1976)  notes  that  1)  storms 
concentrate  along  the  northern  edge  of  the  Gulf  Stream; 

2)  deepening  is  favored  along  the  coast  and  along  the 
northern  edge  of  the  Gulf  Stream;  and  .3)  that  winter  pre¬ 
cipitation  along  the  Atlantic  coast  may  be  modulated 
by  water  temperatures  between  the  Gulf  Stream  and  the  coast. 
The  position  of  the  Gulf  Stream,  the  strength  of  the  baro- 
clinic  zone  along  its  northern  edge  and  the  strength  of 
the  baroclinic  zone  along  the  coast  all  vary  in  time. 
Harrison  et  al.  (1965)  reports  that  coastal  erosion  at  Virginia 
Beach  is  most  highly  correlated  with  water  temperatures. 

This  is  due  to  the  presence  of  storms  offshore  (along  the 
Gulf  Stream)  and  the  high  waves  they  produce  which  in  turn 
erode  the  beach,  mix  the  water  column  and  advect  cold 
waters  southward  along  the  coast.  In  studies  at  the 
University  of  Virginia  we  have  found  that  the  frequency 
maximum  of  cyclogenesis  along  the  coastal  baroclinic 


zone  lags  behind  the  frequency  maximum  of  storms  moving 
along  the  Gulf  Stream  (Resio  and  Hayden,  1975) .  The  intensity 
of  the  coastal  baroclinic  zone  may  be  modulated  by  the 
movement  of  storms  along  the  Gulf  Stream,  and  by  shifts  in 
the  Gulf  Stream  and  eddies  along  its  margins. 

The  winter  geography  of  snowcover  over  North  American 
and  Pacific  Ocean  sea  surface  temperatures  may  also  contribute 
to  the  persistence  of  cyclone  frequency  patterns.  In 
fact,  given  the  multiplicity  of  potential  causes  of  per¬ 
sistence,  it  is  not  surprising  that  several  statistically 
independent  modes  of  cyclone  frequency  persistence  exist. 

Persistence  forecasting  is  well  within  the  tradition 
of  forecast  method  and  technique  and  it  is  generally 
recognized  that  persistence  eventually  fails.  Failure 
(negative  skill)  using  the  University  of  Virginia  Climate 
Prediction  Model  occurs  about  8  percent  of  the  time  for 
winter  forecasts  and  12  percent  for  the  summer.  Failure  may 
occur  because  the  year  forecasted  is  near  the  mean  pattern 
or  because  circulation  patterns  change,  sometimes  abruptly 
(Namias,  1960) .  The  causes  of  the  circulation  change  are 
more  obscure.  Additional  research  on  seasons  with  negative 
forecast  skill  is  required,  before  the  forecast  system 
could  be  applied  with  equal  confidence  in  all  years. 

Preliminary  analyses  indicate  that  forecast  skill  can 
be  improved  by  multiple  regression  models  for  predicting 
eigenvector  weightings.  These  analyses  also  suggest  that 
forecast  skill  will  be  maximum  for  forecasts  slightly 
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longer  than  1  year  in  advance.  New  models  for  these 
improvements  and  for  an  enlarged  geographic  area  of  coverage 
are  being  constructed  and  evaluated  and  will  be  reported  in 
due  course. 
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Figure  7: 


FIGURE  CAPTIONS 

Chart  of  the  study  area.  The  rectangular  inset 
is  2.5°  latitude  by  5.0°  longitude.  There  are 
74  such  rectangular  grid  cells  in  the  study  area. 
The  arrow  represents  a  storm  track  passing 
through  the  grid  cell  shown. 

Mean  Winter  (A)  and  Summer  (B)  cyclone  frequency 
for  the  years  1835-1979.  Winter  includes  the 
months  of  October  thru  March  and  summer  consists 
of  April  through  September.  Units  are  cyclones 
per  2.5°  latitude-5®  longitude  grid  cell. 

Standard  deviations  of  mean  Winter  (A)  and 
Summer  (B)  cyclone  frequency  for  the  years  1885- 
1979. 

Percent  variance  explained  versus  eigenvector 
number  for  the  first  21  eigenvectors  of  Winter 
and  Summer  cyclone  frequency  matrices  and 

for  matrices  of  the  same  dimensions  constructed 
by  randomizing  the  original  winter  and  summer 
matrices . 

The  first  eigenvectors  of  Winter  (A)  and  Summer 
(B)  cyclone  frequencies. 

The  second  eigenvectors  of  Winter  (A)  and  Summer 
(B)  cyclone  frequencies. 

The  third  eigenvectors  of  Winter  (A)  and  Summer 


(B)  cyclone  frequencies. 


Figure  8:  The  fourth  eigenvectors  of  Winter  (A)  and  Summer 


(B)  cyclone  frequencies. 

Figure  9:  Time  variation  of  the  weightings  of  the  first 
four  Winter  (A)  and  Summer  (B)  eigenvectors  of 
cyclone  frequencies.  Eigenvector  number  is 
specified  on  the  ordinates  and  winter  eigen¬ 
vectors  by  the  designation  (A)  and  summer  by 
(B)  . 
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Figure 


10:  Global  percent  skill  scores  for  Winter 
and  Summer  forecasts. 

11:  Histograms  of  frequencies  of  global  percent 
skill  scores  in  5%  class  intervals  for  Winter 
(A)  and  Summer  (B)  forecasts  shown  in  Figure  10. 

12:  Global  percent  skill  scores  versus  global 

quadratic  skill  scores  for  winter  and  summer 
forecasts  combined. 

13:  Summer  and  winter  global  percent  skill  scores 
for  the  1970s  as  calculated  in  Figure  10  'solid 
line)  compared  to  the  same  skill  scores  based 
on  predicting  the  1970s  one  year  at  a  time  (dashed 
line) . 

14:  Global  percent  (solid  line),  Heidke  (long  dash 
line) ,  quadratic  (short  dash  line)  and  deviation 
(dotted  line)  for  the  sequentially  forecasted 
years  of  the  1970s.  Percent  skill  score  is  given 
in  decimal  notation. 

15:  Local  percent  skill  scores  for  Winter  and 


Summer  . 
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Figure  16:  Histograms  of  frequencies  of  local  skill  scores 
in  5%  class  intervals  for  the  Winter  (A)  and 
Summer  (B) .  Forecasts  shown  in  Figure  14. 

Figure  17:  Forecast  for  the  winter  of  1979-1980  and 

verification.  A  shows  the  predicted  cyclone 
frequency  pattern  in  cyclones  per  grid  cell. 

Dashed  lines  indicate  the  predicted  mean  storm 
track  as  evidenced  by  local  maxima  in  the  data 
field.  B  shows  the  actual  occurrences  of  cyclones 
for  the  winter  and  the  mean  storm  track.  C  shows 
the  difference  pattern  between  the  predicted 
and  the  long-term  mean  and  D  the  difference 
between  the  actual  cyclone  occurrence  and  the 
mean.  E  shows  the  difference  between  the  pre¬ 
dicted  and  the  actual.  Solid  circles  indicate 
correct  forecast;  open  circles  are  failures. 


TABLE  1 


Percent  Variance  Explained 


Eigenvector 


Winter 


Summer 


TABLE  2 


Pearson  Correlation 
Coefficients 

.74* 

.52* 

.10 

-.46* 

.74* 

.51* 

.10 

-.46* 


*  =  Significant  at  the  .00001  level 
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TABLE  3 


SUMMARY  OF  ESTIMATES  OF  GLOBAL  FORECAST  SKILL 
For  the  University  of  Virginia  Climate  Prediction  Model 


Winter 

Summer 

1885/86-1978/79 

1886-1979 

PERCENT  SKILL  SCORE 

66% 

66% 

HEIDKE  SKILL  SCORE 

0.34 

0.32 

DEVIATION  SKILL  SCORE 

0.08 

0.10 

QUADRATIC  SKILL  DEVIATION 

0.14 

0.15 

AVERAGE  CYCL*ONES  PER 

GRID  CELL 

5.9 

4.0 

AVERAGE  STANDARD  DEVIATION 
OF  CYCLONES  PER 

GRID  CELL 

2.9 

2.1 

AVERAGE  ABSOLUTE  ERROR 

FOR  FORECAST 

2.30 

1.91 

AVERAGE  ABSOLUTE  ERROR 

FOR  MEAN  AS  FORECAST 

2.57 

2.37 

ROOT  MEAN  SQUARE  ERROR 

FOR  FORECAST 

3.07 

2.55 

FOOT  MEAN  SQUARE  ERROR 

FOR  MEAN  AS  FORECAST 

3.39 

2.86 

TABLE  4 


Summary  of  correlations 
departures  from  fitted 

global  skill  scores 

between  eigenvector  weighting 
third  order  polynomials  (  P  D ) 

(S)  in  the  following  season 

PEARSON 

Correlation 

Regression 

Slope* 

Significance 

Level 

Winter  S  vs  Summer  PD 

E1 

-.14 

-1.30 

.095 

E2 

+  .07 

+  .48 

.246 

E3 

-.05 

-  .49 

.  318 

E4 

+  .03 

+  .40 

.  376 

Summer  S  vs  Winter  PD 

E1 

-.22 

-1.82 

.015 

E2 

-.04 

-  .28 

.343 

E3 

-.04 

-  .40 

.340 

E4 

+  .15 

+  2.31 

.031 

* 


the  units  on  the  regression  slope  are  percentage 
points  of  global  skill  per  unit  eigenvector  weighting 
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ABSTRACT 

To  define  spatial  and  temporal  variations  in  annual  cyclone  frequencies,  principal  components  were 
calculated  from  a  matrix  of  annual  frequencies  for  74  grid  cells  covering  eastern  North  America  and  the 
western  North  Atlantic  and  the  years  1885  to  1978.  The  first  principal  component  contrasts  cyclone 
frequencies  in  continental  versus  marine  areas.  Since  the  early  years  of  this  century,  there  has  been  a 
trend  toward  increased  cyclone  frequency  over  marine  areas  and  a  decline  in  frequencies  over  the  con¬ 
tinent.  This  trend  peaked  in  the  19*0' s.  The  second  principal  component  is  interpreted  as  an  east  coast 
cyclogenesis  function.  Like  the  first  component,  it  exhibits  a  century-long  secular  variation  with  increas¬ 
ing  coastal  cyciogencsis  in  recent  decades  and  a  maximum  in  the  1950's.  The  first  two  components  explain 
45%  of  the  total  variance.  Higher  order  vectors  ( 3rd  and  4th>  explained  1 2%  of  the  variance  and  geographi¬ 
cally  depict  variance  in  (he  Gulf  coast  and  Great  Lakes  regions,  respectively.  Secular  variations  in  the 
weightings  of  (he  third  and  fourth  components  contain  higher  frequency  variations  than  the  first  and 
second  components. 


1.  Introduction 

Damaging  waves  and  storm  surges  along  the  U.S. 
Atlantic  coast  are  largely  due  to  extratropicai 
cyclones.  Along  most  of  this  coast,  shorelines  are 
receding  (Hayden  et  al.,  1979:  Dolan  et  al.,  1979). 
The  erosion  trend  is  attributed  to  the  current  rise  in 
sea  level  (Bruun,  1962).  lower  average  cyclone 
central  pressure  (.Mather  et  al.,  1964),  reduced 
coastal  sand  supply  (Hoyt.  1967),  human  activities 
(Dolan  et  al.,  1973).  and  to  secular  variation  in 
cyclone  frequency,  magnitude  and  duration  (Hay¬ 
den.  1975).  More  recently,  changes  in  the  movement 
of  cyclones  along  the  Atlantic  coast  have  been 
reported  (Resio  and  Hayden,  1975:  Dickson  and 
Namias.  1976).  Changes  in  the  track  of  cyclones 
relative  to  the  coast  give  rise  to  variations  in 
breaker  heights  and  storm  surge  heights  at  the  coast 
(Resio  and  Hayden.  1975). 

The  Resio  and  Hayden  study  is  limited  by  the 
small  reach  of  coast  studied  and  the  Dickson  and 
Namias  study  by  the  short  time  span  analyzed 
(1948-75).  In  this  report,  spatial  and  temporal 
variations  in  cyclone  frequencies  are  examined  for 
the  years  1885  to  1978.  The  study  area  is  eastern 
North  America  and  the  western  North  Atlantic  east 
of  the  100th  meridian  and  south  of  the  50th  parallel 
(Fig.  1). 

2.  Earlier  work  on  cyclone  frequencies 

Pettensen  (1950)  mapped  the  geographical  fre¬ 
quencies  of  cyclones  in  the  Northern  Hemisphere 


using  daily  sea  level  Historical  Weather  Maps  for 
the  period  1899-1939.  Petterssen’s  maps  show  a 
frequency  maximum  along  the  east  coast  of  the 
United  States  as  far  south  as  Florida.  Klein  and 
Winston  (1958)  attribute  the  east  coast  cyclone 
frequency  maximum  to  an  Atlantic  seaboard 
cyclone  track.  Cyclogenesis  associated  with  this 
track  occurs  along  the  tier  of  Gulf  coast  states  or 
along  the  Virginia  Capes.  Reed  (1960)  indicates  in 
his  maps  of  percentage  frequency  of  fronts  that 
the  front-frequency  maximum  also  parallels  the  east 
coast  margin  but  at  some  distance  offshore  during 
the  years  1952-56.  Miller  (1946)  analyzed  cyclo¬ 
genesis  for  the  mid-Atlantic  coastal  region  and 
found  two  distinct  types:  Type  A  which  form 
largely  to  the  north  of  Cape  Hatteras.  North 
Carolina,  and  Type  B  which  form  to  the  south  of 
Cape  Hatteras.  While  seasonal  variations  in  cyclone 
frequency,  front  frequency  and  track  locations  are 
discussed  in  some  detail,  temporal  variations  were 
not  studied. 

Hosier  and  Gamage  (1956)  studied  U.S.  cyclone 
frequencies  in  5°  latitudinal  and  longitudinal  grid 
cells  for  the  50  years  ending  in  1955.  They  found 
no  evidence  in  the  annual  means  of  periodicities, 
trends  or  shifts  in  regions  of  maximum  frequency. 
Reitan  (1974,  1979)  updated  Hosier  and  Gamage's 
work  and  reported  a  general  decline  in  cyclone 
frequencies. 

Recent  studies  of  cyclones  along  the  Atlantic 
coast  have  focused  on  the  generation  of  potentially 
damaging  waves  and  surges.  Hayden  <  1975)  showed 
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an  increase  in  the  frequency  of  storms  that  generate 
deep-water  waves  of  1.3  m  or  greater  from  1942- 
67.  Hayden  concluded  that  the  number,  severity 
and  duration  of  Atlantic  coast  storms  had  increased 
between  1942  and  1974.  Increases  in  cyclone 
severity  are  also  reported  by  Mather  et  ai.  (19641 
and  they  attribute  the  changes  to  decreased  central 
pressure  between  the  1920’s  and  1960’s.  Resio  and 
Hayden  ( 19751  found  that  cyclones  moving  along  the 
U.S.  Gulf  coast  and  the  genesis  of  new  cyclones 
along  the  mid-Atlantic  coast  had  increased  from 
the  1920’s  to  the  1960’s.  The  pre-1920’s  period 
more  closely  resembled  the  1960’s  than  the  years 
following  the  1920’ s.  Resio  and  Hayden  also  showed 
that  the  increased  severity  of  Atlantic  coastal  storms 
resulted  from  a  seaward  displacement  of  mean  storm 
tracks.  They  attributed  this  change  to  increased 
blocking  in  the  high  latitudes.  Dickson  and  Namias 
1 1976)  identified  blocking  in  the  Greenland  area  and 
cooler  than  normal  temperatures  in  the  U.S.  South¬ 
east  as  the  cause  of  increased  frequencies  of 
cyclones  off  the  Atlantic  coast  between  the  1940’s, 
and  the  1950's  and  1960’s. 

3.  The  data 

Cyclone  frequencies  for  each  year  (1885-1978) 
were  tabulated  for  the  74.  2.5°  latitude  by  5.0®  longi¬ 
tude  grid  cells  comprising  the  study  area  shown  in 
Fig.  1,  From  monthly  charts  of  the  "Tracks  of  the 
Centers  of  Cyclones  at  Sea  Level”  published  by  the 
Monthly  Weather  Review  and  in  recent  years  by 


Mariners  Weather  Log.  annual  totals  of  cyclones 
passing  through  each  grid  cell  were  recorded.  Mul¬ 
tiple  entries  of  a  given  storm  in  a  grid  cell  were 
ignored.  Raw  grid-cell  frequencies  were  used  rather 
than  area-normalized  frequencies,  thus  avoiding  a 
latitude-dependent  bias  which  can  occur  as  a  result 
of  area  normalization  (Ballenzweig,  1959).  A  com¬ 
panion  analysis  was  run  on  area-normalized  data. 
The  patterns  of  both  the  long-term  means  and 
standard  deviations  changed  very  slightly  but  the 
calculated  principal  components,  percentages  of 
variance  explained,  and  time  series  of  weightings 
were  unchanged  compared  to  analysis  of  raw  fre¬ 
quency  data.  The  resulting  data  matrix  has  the 
dimensions  of  74  grid-cell  variables  and  94  annual 
cases. 

4.  Long-term  means  and  variances 

Fig.  2  shows  the  mean  annual  frequency  of 
cyclones  for  the  period  1885-1978.  In  general  the 
frequency  of  cyclones  increases  with  latitude  and 
with  proximity  to  the  U.S.  east  coast.  The  axis  of 
the  east  coast  frequency  maximum  is  centered  over 
the  coast  from  Nova  Scotia  to  Georgia.  The  loca¬ 
tion  of  this  frequency  maximum  coincides  with  the 
baroclinic  zone  separating  the  cold  mainland  from 
the  warmer  ocean  waters. 

The  standard  deviation  of  annual  cyclone  fre¬ 
quencies  is  presented  in  Fig.  3.  The  dominant 
feature  shown  by  the  chan  is  the  pronounced 
frequency  maximum  offshore.  The  axis  of  this 


Fig.  I.  Chart  of  the  study  area.  The  rectangular  inset  is  2.3°  latitude  <  5.0° 
longitude.  There  are  74  such  rectangular  grid  cells  in  the  study  area.  The  arrow 
represents  a  storm  track  passing  through  the  grid  cell  shown. 
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Fig.  2.  Mean  annual  cyclone  frequency  for  the  years  1885-1978. 


maximum  is  seaward  of  the  maximum  in  the  mean  margin.  The  standard  deviations  landward  of  the 
field  ( Fig.  2)  and  apparently  rotated  clockwise  margin  of  the  east  coast  are  rather  modest  suggest- 
several  degrees.  The  location  of  the  standard  ing  that  the  axis  of  yearly  storm  occurrences  may 
deviation  maximum  more  closely  approximates  the  expand  eastward  but  in  general  not  over  the 
locus  of  the  Gulf  Stream  than  the  continental  mainland. 


Fig  3.  Standard  deviations  of  the  mean  annual  cyclone  frequency 
for  the  years  1885-1978. 
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Table  1 .  The  percentage  of  the  total  variance 
for  the  Ant  four  eigenvectors. 


Eigenvector 

number 

Percent  variance 
explained 

Cumulative  percent 
variance  explained 

1 

23.0 

28.0 

*> 

17.3 

45.3 

3 

6.6 

51.9 

4 

5.5 

57.4 

5.  Principal  component  analysis  of  annual  frequency 
data 

a.  The  analysis 

Principal  component  analysis  has  successfully 
resolved  the  variance  structure  in  multivariant, 
geophysical  data (Kutzbach.  1967;  Frittseru/.,  1971; 
Resio  and  Hayden.  1975),  and  in  terms  of  least 
square  errors,  this  type  of  analysis  provides  a 
method  for  determining  patterns  in  large  data  fields 
(Lorenz,  1956;  Gilman.  1957;  Kutzbach.  1967).  The 
objective  of  the  analysis  is  to  isolate  characteristic, 
recurrent,  and  independent  modes  of  covariance 
among  variables  into  a  new  set  of  independent 
variables.  Basically,  the  analysis  transforms  a  set  of 
intercorrelated  variables  into  a  new  coordinate 
system  in  which  the  axes  are  linear  combinations 
of  the  original  variables  and  are  mutually  orthogonal. 

To  prevent  those  grid  cells  with  high  mean  cyclone 
frequencies  (high  latitudes)  from  dominating  the 
total  variance  and  consequently  from  dominating  the 
eigenvector  forms,  the  correlation  matrix  was  used 
rather  than  the  covariance  matrix.  The  procedures 
for  calculation  of  the  eigenvectors  follow  those  of 
Kutzbach  (1967)  and  Vincent  et  al.  (1976). 

Principal  component  analysis  provides  a  descrip¬ 
tion  of  the  major  modes  of  variability  in  the  data 
set.  Typically,  each  component  is  identified  with 
some  property  of  the  data  field.  The  analysis  also 
provides  an  index  which  measures  the  importance 
of  each  component  within  each  year.  Finally,  the 
analysis  provides  an  estimate  of  the  total  percent  of 
variance  in  the  data  set  which  can  be  explained 
on  the  basis  of  each  component. 

h.  The  eigenvectors 

The  percentage  of  variance  and  the  cumulative 
percentage  of  variance  explained  by  the  first  four 
eigenvectors  is  given  in  Table  1.  As  indicated  in 
the  table  the  first  four  eigenvectors  account  for 
57.4%  of  the  total  variance.  The  problem  has  thus 
been  reduced  from  a  74-variable  problem  in  which 
the  variables  are  intercorrelated.  to  a  4-variable 
problem  in  which  each  new  variable  is  orthogonal 
and  thus  statistically  independent.  Higher  order 


eigenvectors  explained  smaller  fractions  of  the  total 
variance  and  are  not  dealt  with  here. 

The  eigenvectors  corresponding  to  the  eigen¬ 
values  1  through  4  are  mapped  in  Figs.  4-7.  The 
first  eigenvector  (Fig.  4)  indicates  that  the  dominant 
mode  of  cyclone  frequency  variation  about  the  mean 
is.  with  positive  weightings,  a  relative  abundance  of 
cyclones  over  marine  areas  with  diminished  cyclone 
frequencies  over  the  continent.  Years  with  negative 
weightings  would  be  characterized  by  relatively  few 
cyclones  over  marine  areas  and  above  average 
frequencies  over  the  continent.  The  corridor  of  posi¬ 
tive  values  over  marine  areas  shown  in  Fig.  4  cor¬ 
responds  in  location  to  the  zone  of  increased 
cyclone  frequencies  found  by  Dickson  and  Namias 
( 1976)  and  the  east  and  south  displaced  storm  track 
identified  by  Resio  and  Hayden  (1975).  The  first 
eigenvector  suggests  a  negative  teleconnection  be¬ 
tween  offshore  displacement  of  cyclones  and  the 
frequency  of  Colorado  lows.  Inspection  of  con¬ 
toured  cyclone  frequency  maps  confirms  the  exist¬ 
ence  of  this  negative  teleconnection.  Thus  secular 
variations  in  the  frequency  of  offshore  displace¬ 
ment  of  cyclones  is  out  of  phase  with  the  fre¬ 
quency  of  Colorado  lows. 

The  second  eigenvector  (Fig.  5)  has  positive  val¬ 
ues  over  the  entire  field  and  as  such  indicates  for 
positively-weighted  years,  a  general  increase  in 
cyclone  numbers  and  in  negatively-weighted  years, 
fewer  than  average  cyclones.  The  dominant  feature 
of  the  map  of  eigenvector  no.  2  is  the  maximum 
centered  over  the  coast  of  the  mid-Atlantic.  This 
pattern  is  rather  similar  to  the  charts  of  cyclo- 
genesis  frequency  published  by  Petterssen  (1941  and 
1956)  and  Reitan  (1974).  This  suggests  that  the 
second  principal  component  might  best  be  termed  a 
coastal  cyclogenesis  function. 

The  positive  sense  of  the  third  eigenvector  (Fig. 
6)  contrasts  above-normal  frequencies  of  storms 
moving  eastward  from  the  Great  Plains  up  to  the 
Ohio  Valley  and  exiting  the  coast  between  Cape 
Hatteras  and  Cape  Cod.  with  below  normal  fre¬ 
quencies  of  cyclones  in  the  Gulf  coast  region.  The 
fourth  eigenvector  (Fig.  7)  contrasts  the  frequency 
of  cyclone  occurrence  in  the  Great  Lakes-Saint 
Lawrence  River  region  with  the  frequency  of  cy¬ 
clones  moving  out  of  the  Great  Plains.  In  years 
with  positive  weightings,  storms  are  more  frequent 
than  average  in  southeastern  Canada.  Negatively- 
weighted  years  show  an  increase  in  storm  fre¬ 
quency  in  the  Colorado  storm  track  area. 

The  first  four  eigenvectors  of  annual  cyclone  fre¬ 
quencies  constitute  four  new  orthogonal  axes  which 
account  for  nearly  60%  of  the  variance  in  the 
original  data.  In  general,  weightings  on  these  four 
vectors  for  the  93  individual  years  of  record  varied 
from  -10  to  -10.  The  meanings  of  each  of  the 
four  new  axes  are  summarized  in  Fig.  8. 


eigenvector  of  annual  cyclone  frequencies  Kk,  7.  The  fourth  eigenvector  of  annual  cyclone  frequencies. 
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and  the  -  sign  associated  with  cyclogenesis  for  axis  no.  2  indi¬ 
cates  below  and  above  normal  cyclone  numbers  in  the  region  of 
coastal  cyclogenesis  along  the  U.S.  mid-Atlantic  coast. 


c.  The  time  series 

The  time  series  of  eigenvector  weightings  for  each 
year  (1885-1978)  for  each  of  the  first  four  eigen¬ 
vectors  is  shown  in  Figs.  9- 12.  The  annual  weight¬ 
ings  of  the  first  eigenvector  show  a  broad  quasi- 
sinusoidal  variation  of  period  slightly  greater  than 
100  years.  This  variation  thus  suggests  that  cyclone 
frequencies  declined  over  marine  areas  and  in¬ 
creased  over  the  continent  between  1885  and  1925. 
Since  1925,  cyclone  frequencies  over  marine  areas 
have  increased  and  frequencies  over  the  continent 
have  declined.  With  the  limited  time  series  length 
it  is  not  possible  to  establish  whether  the  variation 
is  indeed  periodic. 

Annual  weightings  of  the  second  eigenvector 
(Fig.  10)  also  exhibit  a  century-long  scale  of  varia¬ 
tion  from  negative  values  at  the  beginning  of  the 
time  series  to  positive  values  after  1910  and  appar¬ 
ently  back  to  negative  values  around  1960.  If  the 
interpretation  of  eigenvector  2  as  representing 
coastal  cyclogenesis  is  correct,  then  there  should 
have  been  an  increase  in  coastal  cyciogenesis  during 
the  first  four  decades  of  the  20th  century.  This 
conclusion  is  consistent  with  the  work  of  Hosier  and 
Gamage  ( 1956)  and  Reitan  ( 1979).  The  secular  varia¬ 
tion  shown  in  Fig.  9  is  remarkably  similar  to  that 
shown  by  Reitan  (1979)  in  annual  cyclone  totals  (his 
Fig.  3).  Reitan  (1979)  notes  that  about  30%  of  the 
decline  in  cyclone  frequency  since  1950  is  due  to 
decreased  cyciogenesis.  The  early  1940’ s  minimum 
of  cyclone  totals  reported  by  Reitan  (1979)  is 
clearly  evident  in  the  time  series  weighting  of  the 
second  principal  component  (Fig.  9). 

The  secular  variation  evident  in  Figs.  9  and  10 
raises  questions  about  the  detectability  of  cyclones 
off  the  U.S.  Atlantic  coast  in  the  early  portion  of  the 
record.  Examination  of  year-by-year  cyclone  fre¬ 


quency  charts  shows  a  scarcity  of  cyclones  over 
the  continent  of  the  same  order  in  recent  decades 
as  over  the  oceanic  areas  earlier  in  the  century. 
Resio  and  Hayden  (1975)  used  1899-1970  hemi¬ 
spheric  pressure  data  to  reconstruct  mid-Atlantic 
coast  cyclone  frequencies.  This  synthetic  record  for 
cyclones  migrating  from  inland  and  forming  along 
the  coast  shows  a  minimum  in  the  1920's  and  a 
maximum  in  the  1960’s  similar  to  those  shown  in 
Figs.  9  and  10. 

To  determine  whether  an  observational  insuf¬ 
ficiency  in  marine  areas  in  the  early  years  of  record 
would  account  for  the  secular  trend  shown  in  Fig.  9 
a  reduced  data  matrix  was  analyzed.  All  grid  ceils 
south  of  27.5°N  and  east  of  70°W  were  removed 
and  the  principal  components  were  calculated. 
There  were  essentially  no  differences  in  the  result¬ 
ing  components  or  time  series  of  weightings.  The 
percentage  of  variance  explained  by  the  first  four 
components  increased  2%.  The  first  eigenvector  for 
the  reduced  data  matrix  is  shown  in  Fig.  13.  The 
similarity  to  Fig.  4  is  unmistakable.  Clearly  the 
shape  of  the  eigenvector  is  not  changed  by  chang¬ 
ing  the  size  of  the  data  matrix  or  the  removal  of 
marine  grid  cells  from  the  data.  There  seems  to  be  no 
systematic  observational  insufficiency  associated 
with  the*marine  areas  of  the  data  field.  The 
highest  weightings  in  the  first  vector  (Fig.  4)  are  not 
found  in  the  middle  of  the  grid  as  suggested  by 
Buell  (1975).  Changing  the  grid  size  did  not  alter 
the  pattern.  The  secular  trends  then  must  be  con¬ 
sidered  real  and  not  observational  or  mathematical 
artifacts. 

The  time  series  of  the  third  and  fourth  eigen¬ 
vectors  (Fig.  11  and  Fig.  12)  show  excursions 
from  positive  to  negative  weightings  over  shorter 
periods  of  years.  The  third  eigenvector  runs  from 
negative  values  to  positive  values  between  1900  and 
1925,  1925  and  1950  and  between  1955  and  1978. 
Spectral  analysis  shows  some  concentration  of 
spectral  power  at  about  23  years.  The  fourth 
eigenvector  shows  a  broad  decline  from  positive 
values  in  the  early  part  of  the  record  to  low 
values  in  the  late  1940's.  Since  the  1940’s  weight¬ 
ings  have  again  increased. 

The  original  data  matrix  was  divided  into  a  winter 
(October-  March)  and  a  summer  <  April- September) 
half.  Principal  components  were  calculated.  This 
analysis  was  done  to  insure  that  the  annual  cyclone 
frequency  principal  components  were  not  meaning¬ 
less  summations  of  winter  and  summer  conditions. 
The  forms  of  each  set  of  principal  components 
were  very  similar  to  those  shown  in  Figs.  4  through 
7  and  the  time  series  weightings  were  remarkably 
similar  to  those  shown  in  Figs.  9  through  12.  It 
appears  that  the  main  time  scale  of  variation  is 
extra-  rather  than  inter-annual. 


weightings  of  the  third  eigenvector  (fc'4).  h'Ki  12.  Time  variii'ion  of  the  annual  weightings  of  the  fourth  eigenvector  (£,). 
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Ftc.  13.  The  first  eigenvector  of  annual  cyclone  frequencies  for  a  reduced  data 
matrix.  The  latitude  limits  of  the  data  held  are  27.5°N  and  47.5'N.  The  longitude 
limits  are  70*W  and  lOtTW. 


6.  Discussion  and  conclusion 

The  variance  in  cyclone  frequency  is  maximum 
off  the  U  S.  Atlantic  coast.  Much  of  this  variance 
is  due  to  the  eastward  displacement  of  storm  tracks 
that  parallel  the  coast  and  to  cyclogenesis  along  the 
mid-Atlantic  portion  of  the  east  coast.  Both  of  these 
two  contributions  to  the  total  variance  exhibit  cen¬ 
tury-long  secular  variation.  These  changes  are  un¬ 
doubtedly  due  to  changes  in  the  intensity  of  the 
east  coast  baroclinic  zone  and  to  shifts  in  the  North 
American  long- wave  location  associated  with  block¬ 
ing  in  the  high  latitudes  as  suggested  by  Resio  and 
Hayden  (1975)  and  Dickson  and  Namias  (1976). 

It  thus  appears  that  the  increased  storminess  of 
the  U.S.  east  coast  since  the  early  1940’s  (Mather 
et  ai,  1964;  Hayden.  1975:  Resio  and  Hayden,  1 97 5 ■_ 
Dickson  and  Namias,  1976)  is  part  of  a  secular  varia¬ 
tion  of  a  longer  time  scale.  While  the  motivation  for 
this  study  was  the  apparent  increase  in  east  coast 
damaging  cyclones  during  this  century,  the  results 
of  the  study  have  wider  significance.  In  general,  we 
may  conclude  that  there  have  been  several  secular 
variations  in  cyclone  frequency  patterns:  1)  declines 
in  the  frequency  of  Colorado  lows  have  been  as¬ 
sociated  with  increases  in  cyclones  off  the  Atlantic 
coast:  2)  a  general  decline  in  cyclone  frequencies 
is  accompanied  by  a  decline  in  cyclogenesis  along 
the  mid- Atlantic  coast. 
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ABSTRACT 

Meteorological  data  extracted  from  maps  and  charts  is  frequently  normalized  for  area  because  grid  cell 
size  changes  with  latitude.  Raw  and  area- normalized  cyclone  frequency  data  were  compared.  .Area 
normalization  deforms  mean  frequency  patterns.  The  patterns  of  standard  deviations  about  the  mean 
are  not  deformed  in  the  same  way.  Area  normalization  does  not  change  calculated  principal  components 
or  their  time  histories  of  weightings. 


1.  Introduction 

Bailenzweig  (1959)  detailed  the  advantages  of 
equal-area  grids  in  data  extraction  and  in  mapping 
the  geographical  distribution  of  physical  quantities. 
He  also  listed  the  disadvantages  of  using  latitude - 
longitude  grid  ceils  and  of  the  subsequent  use  of 
latitude-dependent  correction  to  achieve  an  area 
normalization.  Ballenzweig’s  "practical”  equal- 
area  grid  has  not  been  widely  used  in  data  extrac¬ 
tion  from  maps  in  spite  of  its  merits.  This  has  oc¬ 
curred.  ironically,  due  to  practical  considerations 
arising  from  the  diversity  of  scale  and  projection 
of  historical  data  in  map  form.  For  detailed  collec¬ 
tion  from  maps,  preparation  of  new  equal-area  grids 
each  time  scale  projection  changes  is  a  major  dis¬ 
advantage.  In  almost  all  cases  data  are  extracted 
from  published  maps  according  to  grid  cells  defined 
by  lines  of  latitude  and  longitude. 

Judging  from  the  recent  literature  (O'Connor. 
1964;  Taljaard.  1967;  Reitan.  1974.  1979),  latitude- 
dependent  area  normalization  has  become  the 
method  of  choice  even  though  Bailenzweig  discour¬ 
ages  its  use  in  favor  of  raw  frequencies  in  n  degree 
squares.  South  of  the  reference  latitude,  area- 
normalization  adjusts  the  frequency  downward,  to 
make  it  represent  a  box  smaller  than  the  one  in 
which  the  original  tabulations  were  made.  The  fre¬ 
quency  of  cyclones  is  thereby  artificially  decreased. 
North  of  the  reference  latitude  a  reverse  problem 
exists,  and  cyclone  frequencies  are  artificially  in¬ 
creased.  In  a  recent  study  on  cyclone  frequencies 
(Hayden,  1980)  we  examined  the  issue  anew.  We, 
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like  most  investigators,  did  not  construct  equal-area 
grid  cells  for  the  data  extraction  process.  This 
decision  was  one  of  convenience.  The  diversity  on 
map  formats  for  published  tracks  of  cyclone  centers 
between  1885  and  1979  weighed  against  constructing 
equal-area  data  collection  grids.  We  have  com¬ 
pared  results  of  an  analysis  using  area- normalized 
data  with  results  using  raw  frequency  totals,  and 
found  no  significant  differences  between  the  two. . 


2.  The  test  data 

Cyclone  frequencies  for  summers  i  April  through 
September)  for  each  year  (1885-1979)  were  tabu¬ 
lated  for  the  74  2.5°  latitude  by  5.0°  longitude  grid 
cells  comprising  the  study  area  shown  in  Fig.  1. 
From  monthly  charts  of  the  "Tracks  of  the  Centers 
of  Cyclones  at  Sea  Level"  published  by  the  Monthly 
Weather  Review  and  in  recent  years  by  Mariners 
Weather  Log,  summer  season  totals  of  cyclones 
passing  through  each  grid  cell  were  recorded.  Mul¬ 
tiple  entries  into  grid  cells  were  ignored.  Two  test 
data  matrices  were  prepared.  One  matrix  consisted 
of  the  raw  grid  cell  frequency  data.  The  second 
matrix  was  area  normalized  for  an  area  103  km- 
using  latitude-dependent,  degree-to-length  conver¬ 
sions  published  in  Table  163  of  the  Smithsonian 
Meteorological  Tables.  In  this  study  the  reference 
latitude  was  north  of  ail  grid  cells  and  thus  fre¬ 
quencies  are  reduced  at  all  locations.  All  grid  ceils 
in  the  study  area  are  slightly  larger  than  10*  km* 
with  the  largest  departures  in  the  lowest  latitudes. 
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Fig.  1.  Chan  of  the  study  area.  The  rectangular  inset  is  2.5°  latitude  by  5.0* 
longitude.  There  are  74  such  rectangular  grid  cells  in  the  study  area.  The  arrow 
represents  a  storm  traclc  passing  through  the  grid  cell  shown. 


'Long-term  means  and  standard-deviations  of  each 
matrix  were  calculated. 

3.  The  means 

Fig.  2  shows  the  means  of  summer  cyclone  fre¬ 
quencies  for  the  raw  frequency  data,  the  area-nor¬ 
malized  frequencies  and  the  difference  between  the 
two  mean  charts.  General  inspection  of  the  charts 
reveals  little  difference  in  the  two  means  with  the 
exception  of  slightly  lower  values  in  the  area-nor¬ 
malized  chart.  The  difference  chan  (Fig.  2c)  shows 
that  the  deformation  of  the  raw  frequency  caused  by 
area  normalization  does  not  vary  smoothly  with  lati¬ 
tude.  The  greatest  difference  is  near  latitude  40° 
where  the  raw  frequency  means  are  generally  large. 

4.  The  standard  deviations 

Fig.  3  shows  the  standard  deviations  associated 
with  the  mean  summer  cyclone  frequencies  shown 
in  Fig.  2.  The  general  pattern  of  the  standard  devia¬ 
tions  of  the  raw  and  area-normalized  frequency  data 
are  similar.  As  expected  the  standard  deviations  are 
slightly  smaller  everywhere.  The  difference  between 
the  two  standard  deviation  chans  again  shows  that 
the  largest  differences  are  near  40°  latitude  with  the 
greatest  reduction  in  the  magnitude  of  the  standard 
deviation  occurring  where  the  mean  values  tend  to 
be  highest.  Again  a  simple  latitudinal  adjustment 
of  the  data  results  in  a  relatively  complex  deforma¬ 
tion  of  the  raw  frequency  data.  The  difference  pat¬ 


terns  for  the  two  means  and  two  standard  deviations 
are  not  identical. 

5.  Discussion 

The  pattern  deformation  caused  by  an  area  nor¬ 
malization  of  the  raw  frequency  data  is  in  part  de¬ 
pendent  on  the  reference  area  of  normalization.  In 
this  study  10s  km2  was  selected  because  it  was  the 
best  round  number  area  closest  to  the  actual  size  of 
the  latitude-longitude  grid  cells  used  in  the  original 
data  collection.  This  choice  is  common:  Reitan 
(1974)  used  5.5  x  10s  km2;  Petterssen  (1950)  chose 
10*  km2;  Keegan  (1958)  used  10s  mi2.  Klein  (1958). 
Taljaard  (1967),  and  O’Connor  (1964)  used  a  refer¬ 
ence  latitude-longitude  grid  cell  at  a  specified  lati¬ 
tude  as  their  standard  for  area  normalization. 

Generally  only  the  average  frequency  charts  are 
presented  and  discussed.  The  deformation  of  mean 
patterns  caused  by  area  normalization  would  not 
seriously  afreet  their  interpretations.  Equal-area 
adjustments  do  change  frequency  totals.  This  prob¬ 
lem  becomes  serious  if  the  area-normalized  fre¬ 
quencies  are  plotted  on  unequal  area  map  projec¬ 
tions  as  is  usually  the  case. 

We  have  run  tests  and  have  found  that  area  nor¬ 
malization  results  in  no  alteration  of  the  principal 
components  eigenvectors  or  the  time  histories  of 
eigenvector  weightings  published  earlier  (Hayden. 
1980).  This  is  because  the  normalization  of  the  data 
inherent  in  the  formation  of  the  correlation  matrix, 
from  which  the  eigenvectors  were  extracted,  cancels 


ence  between  area  normalized  and  raw  lrci|ueucics>. 
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the  effects  of  any  latitude-dependent  scaling  of  the 
variables  involved. 

It  is  concluded  that  equal-area  grids  for  data  ex¬ 
traction  are  desirable,  but  unfortunately  they  are 
impractical.  Latitude-dependent  area  adjustments, 
made  after  the  data  extraction  is  completed,  intro¬ 
duce  a  latitude-dependent  bias  in  the  frequency  pat¬ 
terns.  Such  latitude-dependent  adjustments  have 
been  made  by  many  investigators,  but  probably  do 
not  affect  their  results  enough  to  alter  their  con¬ 
clusions.  However,  it  is  recommended  that  area- 
normalization  be  avoided  in  future  studies,  because 
it  introduces  a  systematic  bias  in  the  analysis. 
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